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NOTICE 


The Council has pleasure in announcing publication, early in 1961, of an 
extra volume of the Geophysical Journal (Volume 4), in honour of the seventieth 


birthday of 
Sir Harold Jeffreys 


The Geophysical Editors have obtained the co-operation of many of Sir 
Harold’s distinguished colleagues and pupils, and believe that the contents of 
this special volume will represent the most advanced work in a number of 
fields in which Sir Harold has long been interested. 
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The complete volume, suitably bound in cloth boards, will also be available 


as a special publication. 
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Orders, with remittance, should be sent to: 
Assistant Secretary, Royal Astronomical Society, Burlington House, London, W.1 


SPECIAL NOTICE TO FELLOWS 


Fellows of the Society may obtain this bound volume of the Geophysical 
Journal at a specially reduced price of £1 tos. od. (US$ 4.50). It is particularly 
requested that Fellows wishing to receive a copy should send their orders, with 
remittance, as soon as possible to: 


Assistant Secretary, Royal Astronomical Society, Burlington House, London, W.1 
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Presentation of Papers at Meeting = 


At some meetings of the Society the background and conclusions of selected papers 
tod In order to assist the Secretaries in the selection of = = 
papers for such meetings, authors are asked to let the Society know, when submitting = = - 
papers, whether they would be willing to give an account of their paper, if requested. 

The attention of authors resident ab is drawn to the fact that the Society welcomes 
information about their work. The Secretaries would be happy to consider havingsuch = 


1. General,—It is the aim of the Society to be of the greatest possible service in = 
disseminating astronomical results and ideas to the scientific community with the utmost 
possible speed. Contributors are accordingly urged to give the most careful considera- 
tion to the presentation of their work, for attention to detail will assuredly result ina ©. 

submitted for publication in Monthly Notices; experience has shown that hreansnely OT ae 
comments of referees have enabled authors to improve the presentation of their work © 
and so increase its scientific value. 4 

2. Communication.—Papers must be communicated to the Society by a Fellow. 

They should be accompanied by a summary at the beginning of the paper conveying 
briefly the content of the paper, and drawing attention to important new information 
and to the main conclusions. The summary should be intelligible in itself, without  — 
reference to the paper, to a reader with some knowledge of the subject; it should not _- 
normally exceed 200 words in length. Authors are requested to submit MSS. in _— 
duplicate. These should be typed using double spacing and leaving a i 
less than one inch on the left-hand side. Corrections to the MSS. s 

the text and not in the margin. i isi 


4- Notation.—For technical astronomical terms, authors should conform closely to 
the recommendations of Commission 3 of the International Astronomical Union (Trans. _ 
I.A.U.; Vol. VI, p. 345, 1938). Council has decided to the LA.U. 3- 


abbreviations for ns where contraction is desirable (Vol. IV, p. 221, 1932), 
In general matters, authors should follow the recommendations in Symbols, ens and | 
Abbreviations (London: Royal Society, 1951) except where these conflict with I.A.U. 
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by the author they are destroyed. 
ne 3. Presentation.—Authors are allowed considerable latitude, but they are requested sii “a aes 
to follow the general style and arrangement of Monthly Notices. Referencestoliterature = -_-s- & 
should be given either in the traditional form of a numbered list at the end of the paper, = _ ; 
or as prescribed in Notes om the Preparation of Papers to be Communicated to the Royal == (sd 
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should Blocks are 10 years; unless the author 
requires before the end o destroyed. Rough copies = “st 
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7. Proofs.—Authors arc liable for costs of alteration exceeding § per cent of = == 
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. maximum conciseness and simplification of symbols and equations consistent with clarity. =| 
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BIBLIOGRAPHICAL NOTE 


The series Occasional Notes of the Royal Astronomical Society, of which the first 


number appeared in 1938 June, terminates with Number 21, ending Volume 3. 


Papers of the same nature as those which have in the past been published in Occasional 
Notes will in future be incorporated in the Quarterly Journal of the Royal Astronomical 
Society, first issued in 1960 September. 
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ROYAL ASTRONOMICAL SOCIETY 
Vol. 121 No. 5 


THE SOLAR RED-SHIFT 


L. A. Higgs 
(Communicated by the Director, University Observatory, Oxford) 
(Received 1960 July 27) 


Summary 

The centre-to-limb increase in wave-length of three medium strength 
iron lines in the solar spectrum has been determined, using micrometer 
measurements of grating spectra. ‘The observations were made at solar 
latitudes of approximately +55°, and at radial distances from the centre 
between 85 and 100 per cent of the radius of the disk. The limb-effect 
curve obtained, after corrections for scattering of light, is in good agreement 
with the curve obtained recently at Oxford (5) for solar equatorial regions. 

The absolute wave-lengths in the solar spectrum of the lines used in this 
investigation and the vacuum arc wave-lengths of the three iron lines were 
determined using the method of circular channels. The centre red-shifts 
of the three lines, together with the limb-effect curve, give a red-shift at 
the extreme limb substantially exceeding the shift required by the principle 
of equivalence. 

It is shown that the line profiles of the solar lines in the grating spectra 
are affected by the apparatus function of the spectrograph. Allowing for 
this, the solar lines at the centre appeared to be symmetric but an 
unexplainéd asymmetry at the extreme limb was discovered. This indicates 
that wave-length measurements alone are insufficient for a determination of 
the nature of limb effect. 


Introduction 


It is an established fact that the wave-lengths of the Fraunhofer lines in the 
solar spectrum are dependent upon the point of observation on the solar disk. 
In general, the wave-length of a line increases as the point of observation moves 
toward the limb. This difference of wave-length between the limb and the centre 
of the disk has received the general name ‘‘limb effect’’. In addition to this 
effect, most Fraunhofer lines have longer wave-lengths than do the corresponding 
transitions observed in a terrestrial laboratory. This effect will be referred to as 
the ‘‘ red-shift’ of the spectral lines. Thus the red-shift of a line at a given disk 
position is given by the red-shift of the line at the centre plus the limb effect. 

The principle of equivalence predicts a wave-length shift of all solar lines 
with respect to corresponding terrestrial lines observed under similar physical 
conditions by an amount given by 


AA 
2°12 x 107%, 


This shift, referred to as the ‘‘ gravity shift’’, is, however, independent of the 
point of observation on the Sun. Thus other physical processes must be invoked 
to describe the observed variations in the red-shifts at the centre of the disk and 
for the very existence of limb effect. 
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Summary 

The centre-to-limb increase in wave-length of three medium strength 
iron lines in the solar spectrum has been determined, using micrometer 
measurements of grating spectra. ‘The observations were made at solar 
latitudes of approximately +55°, and at radial distances from the centre 
between 85 and 100 per cent of the radius of the disk. The limb-effect 
curve obtained, after corrections for scattering of light, is in good agreement 
with the curve obtained recently at Oxford (5) for solar equatorial regions. 

The absolute wave-lengths in the solar spectrum of the lines used in this 
investigation and the vacuum arc wave-lengths of the three iron lines were 
determined using the method of circular channels. The centre red-shifts 
of the three lines, together with the limb-effect curve, give a red-shift at 
the extreme limb substantially exceeding the shift required by the principle 
of equivalence. 

It is shown that the line profiles of the solar lines in the grating spectra 
are affected by the apparatus function of the spectrograph. Allowing for 
this, the solar lines at the centre appeared to be symmetric but an 
unexplained asymmetry at the extreme limb was discovered. This indicates 
that wave-length measurements alone are insufficient for a determination of 
the nature of limb effect. 


Introduction 


It is an established fact that the wave-lengths of the Fraunhofer lines in the 
solar spectrum are dependent upon the point of observation on the solar disk. 
In general, the wave-length of a line increases as the point of observation moves 
toward the limb. This difference of wave-length between the limb and the centre 
of the disk has received the general name ‘‘limb effect’’. In addition to this 
effect, most Fraunhofer lines have longer wave-lengths than do the corresponding 
transitions observed in a terrestrial laboratory. This effect will be referred to as 
the ‘‘ red-shift ’’ of the spectral lines. Thus the red-shift of a line at a given disk 
position is given by the red-shift of the line at the centre plus the limb effect. 

The principle of equivalence predicts a wave-length shift of all solar lines 
with respect to corresponding terrestrial lines observed under similar physical 
conditions by an amount given by 


= = 2°12 x 10-8, 


This shift, referred to as the ‘‘ gravity shift’’, is, however, independent of the 
point of observation on the Sun. Thus other physical processes must be invoked 
to describe the observed variations in the red-shifts at the centre of the disk and 
for the very existence of limb effect. 
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The early history of red-shift and limb-effect determinations may be found 
in the literature (14, 12, 7, 20, 13, 11). In 1948, M. G. Adam (1) published 
the results of an investigation of the wave-lengths of several solar lines at various 
positions along a polar disk radius. It was found that the limb effect reached 
an appreciable magnitude only along the last 10 per cent of a solar radius, and 
indicated the necessity of more detailed observations near the limb. Recently, 
in 1959, such a redetermination of solar limb effect was published by M. G. Adam 
(5). This gave the results of a study of the centre-limb variation of the 
wave-lengths of three medium strength iron lines in solar equatorial regions. 
The wave-length shifts were determined from micrometer measurements of 
spectra obtained using the 32 cm image of the larger Oxford telescope and its 
associated grating spectrograph. ‘These indicated that the red-shift observed at 
the limb for these lines much exceeded the gravity-shift. Recent observations by 
Schréter (19) of limb effect agree well with the Oxford measurements in the 
common range of observation. 

Because of the great theoretical importance of the shift at the limb, a further 
determination of the limb-etfect curve for the same three solar lines as used 
in (§) has been carried out, and forms the subject of this paper. Although no 
great differences exist between the 1948 and 1959 limb-effect curves, referring 
respectively to polar and equatorial regions, the present determination was 
made at intermediate solar latitudes (~ + 55°) with a view to the detection of 
a latitude effect if such exists. 

Part I presents the micrometer determination of the centre-limb variation 
of the wave-lengths of the three solar lines following the procedure adopted by 
M.G. Adam for the 1959 determination. In Part II, the determination of the 
absolute wave-lengths at the centre of the disk of the lines used in the investigation 
is given. Part III deals with an apparent discrepancy arising from the measure- 
ments in Parts I and II and presents some preliminary measurements of the 
profiles of the three solar lines at the centre and limb. 


Part | 

Limb effect from micrometer measurement 

1.1. Method.—The same method of determining the centre-limb change in 
wave-length was adopted as was used in (§). The method is a differential one 
in which the positions of the solar lines on a spectral plate are compared with 
those of a set of standard lines, in this case terrestrial oxygen lines, whose 
wave-lengths are assumed constant. Since observations were confined to the 
part of the disk near the limb and the centre of the disk, such a method requires 
two essentially different sets of measurements. First, a set of standard positions 
of the solar lines at the centre of the disk with respect to the reference lines is 
required, and secondly, a set of measurements of the positions of the solar lines 
at various disk positions with respect to those of the standard lines is needed. 
Then the difference between the two gives the limb effect at various disk positions, 
when corrections have been made for the effects of the rotation of the Sun and 
motion of the observer. These shifts can then be converted to a limb effect 
expressed either in wave-length units or, as is customary, in equivalent Doppler 
velocities. For such a procedure to be valid, the dispersion used in the two sets 
of measurements must be identical, otherwise corrections must be applied, all 
measurements being reduced to one standard dispersion. 
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We shall use the following terminology for disk positions. A given radial 
position will be specified by 6, where @ is the angle between the line of sight 
and a true solar radius at the point of observation. Alternatively, a position is 
specified by p/p, where pp is the radius of the solar disk, and p is the radius to the 
point of observation. To a high degree of approximation, these are related by 

sin = p/py. 

1.2. Solar spectra.—The spectral plates used in the determination were some 
of a series obtained in 1957 (5), others of which were used in the 1959 limb-effect 
study. They were taken using the 35 m solar telescope with its associated 
spectrograph, the Bausch and Lomb grating giving a dispersion of 1-4mm/A 
in the first order. The plates used were all Ilford panchromatic type R50 and 
were developed in Parkhurst two-solution developer. 

Ten plates taken at the centre of the disk were used for the determination of 
standard positions and six plates of the region near the limb were used. For 
the limb plates, the solar image was guided so that the spectrograph slit was 
radial, the limb falling upon the slit. These plates were divided equally between 
the east and west limbs. The sixteen plates used are listed in Table I. 


TABLE | 
Details of measured plates 
Centre plates Limb plates ' 
Date U.T. 1957 Date U.T. 1957 Limb Solar latitude 


May 21°3812 May 24°3929 —50°0 to —59'0 
22°6342 
27°4411 24°3979 — 48-7 to —57°8 
28-4674 
31°4592 27°3944 48:1 to 60°0 
June 12°4850 
13°5754 June 16-4084 — §0°6 to —63°7 
13°6131 
14°4505 16°4141 61°8 
14°6085 
16°4233 East 48-9 to 
The spectral lines are the same as those used by M. G. Adam and are listed 
in Table II, with their I.A.U. wave-lengths. There is no I1.A.U. wave-length 
for the third oxygen line and the value given is an estimated one (§). Table II 
also gives their Rowland intensities, and for the solar lines, the excitation potentials 
of the lower levels and the transitions involved. The oxygen lines are associated 
with the « band of atmospheric oxygen. 


Tasie II 
Measured lines 

1.A.U. Rowland Excitation 
wave-length intensity potential 
6295°960 A 
6297°799 2°21 volts 
6299°228 
6301508 3°64 
(6302-000) 
6302°499 3°67 
6302°764 


Transition 
a®P-y5D® 
25 
31° 


422 L. A. Higgs Vol. 12 


The early history of red-shift and limb-effect determinations may be found 
in the literature (14, 12, 7, 20, 13, 11). In 1948, M. G. Adam (1) published 
the results of an investigation of the wave-lengths of several solar lines at various 
positions along a polar disk radius. It was found that the limb effect reached 
an appreciable magnitude only along the last 10 per cent of a solar radius, and 
indicated the necessity of more detailed observations near the limb. Recently, 
in 1959, such a redetermination of solar limb effect was published by M. G. Adam 
(5). This gave the results of a study of the centre-limb variation of the 
wave-lengths of three medium strength iron lines in solar equatorial regions. 
The wave-length shifts were determined from micrometer measurements of 
spectra obtained using the 32 cm image of the larger Oxford telescope and its 
associated grating spectrograph. ‘These indicated that the red-shift observed at 
the limb for these lines much exceeded the gravity-shift. Recent observations by 
Schréter (19) of limb effect agree well with the Oxford measurements in the 
common range of observation. 

Because of the great theoretical importance of the shift at the limb, a further 
determination of the limb-effect curve for the same three solar lines as used 
in (5) has been carried out, and forms the subject of this paper. Although no 
great differences exist between the 1948 and 1959 limb-effect curves, referring 
respectively to polar and equatorial regions, the present determination was 
made at intermediate solar latitudes (~ + 55°) with a view to the detection of 
a latitude effect if such exists. 

Part I presents the micrometer determination of the centre-limb variation 
of the wave-lengths of the three solar lines following the procedure adopted by 
M.G. Adam for the 1959 determination. In Part II, the determination of the 
absolute wave-lengths at the centre of the disk of the lines used in the investigation 
is given. Part III deals with an apparent discrepancy arising from the measure- 
ments in Parts I and II and presents some preliminary measurements of the 
profiles of the three solar lines at the centre and limb. 


Part | 
Limb effect from micrometer measurement 


I.1. Method.—The same method of determining the centre-limb change in 
wave-length was adopted as was used in (§). ‘The method is a differential one 
in which the positions of the solar lines on a spectral plate are compared with 
those of a set of standard lines, in this case terrestrial oxygen lines, whose 
wave-lengths are assumed constant. Since observations were confined to the 
part of the disk near the limb and the centre of the disk, such a method requires 
two essentially different sets of measurements. First, a set of standard positions 
of the solar lines at the centre of the disk with respect to the reference lines is 
required, and secondly, a set of measurements of the positions of the solar lines 
at various disk positions with respect to those of the standard lines is needed. 
Then the difference between the two gives the limb effect at various disk positions, 
when corrections have been made for the effects of the rotation of the Sun and 
motion of the observer. These shifts can then be converted to a limb effect 
expressed either in wave-length units or, as is customary, in equivalent Doppler 
velocities. For such a procedure to be valid, the dispersion used in the two sets 
of measurements must be identical, otherwise corrections must be applied, all 
measurements being reduced to one standard dispersion. 
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We shall use the following terminology for disk positions. A given radial 
position will be specified by 6, where @ is the angle between the line of sight 
and a true solar radius at the point of observation. Alternatively, a position is 
specified by p/p, where pp is the radius of the solar disk, and p is the radius to the 
point of observation. To a high degree of approximation, these are related by 


sin = p/pp. 


1.2. Solar spectra.—The spectral plates used in the determination were some 
of a series obtained in 1957 (5), others of which were used in the 1959 limb-effect 
study. They were taken using the 35m solar telescope with its associated 
spectrograph, the Bausch and Lomb grating giving a dispersion of 1-4mm/A 
in the first order. The plates used were all Ilford panchromatic type R50 and 
were developed in Parkhurst two-solution developer. 

Ten plates taken at the centre of the disk were used for the determination of 
standard positions and six plates of the region near the limb were used. For 
the limb plates, the solar image was guided so that the spectrograph slit was 
radial, the limb falling upon the slit. These plates were divided equally between 
the east and west limbs. The sixteen plates used are listed in Table I. 


TaBLe 
Details of measured plates 
Centre plates Limb plates 
Date U.T. 1957 Date U.T. 1957 Limb Solar latitude 
May 21°3812 May 24°3929 —50°'0 to —59'0 
22°6342 
27°4411 24°3979 — 48-7 to —57°8 
28-4674 
31°4592 27°3944 48-1 to 600 
June 12°4850 
13°5754 June 16-4084 —§0°6 to —63°7 
13°6131 
14°4505 16-4141 61°8 
14°6085 
16°4233 East 48-9 58°5 
The spectral lines are the same as those used by M. G. Adam and are listed 
in Table II, with their I.A.U. wave-lengths. There is no I.A.U. wave-length 
for the third oxygen line and the value given is an estimated one (§). Table II 
also gives their Rowland intensities, and for the solar lines, the excitation potentials 
of the lower levels and the transitions involved. The oxygen lines are associated 
with the « band of atmospheric oxygen. 


TABLE II 
Measured lines 
1L.A.U. Rowland Excitation 
wave-length intensity potential 
6295°960 A . 3 


6297°799 5 2°21 volts 
6299°228 


6301°508 3°64 
(6302-000) 

6302°499 3°67 
6302°764— 


31* 
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1.3. Centre measurements.—The spectra were measured on a Hilger photo- 
measuring micrometer fitted with a transverse micrometer motion. A diaphragm 
in the focal plane of the eyepiece limited the field of view to a strip of spectrum 
0°145 mm wide. ‘This was positioned parallel to the plate dispersion so that 
only a small section of each spectral line was visible. The plates were measured 
in the forward and reverse directions at the slit height corresponding to the centre 
of the disk and average line positions were found. These were referred to the first 
oxygen line after corrections for the observer’s motion had been applied to the solar 
lines, the dispersions for the three lines being 33-99, 33°96 and 33-95 km/sec/mm. 
The average positions of the seven lines were then found from the ten plates. 
These were adopted as the standard centre positions and the oxygen positions 
were used to indicate the standard dispersion. These standard centre positions 
are listed in Table VI, Part III, where they are compared with the standard 
positions used by M. G. Adam (5). 

1.4. Limb measurements.—Each limb plate was measured at about forty slit 
heights, each corresponding to a different radial position (different p/p)). The 
positions of measurement were so chosen that they fell at approximately equal 
intervals in cos@. In order to compute the value of p/p, for each ‘‘ spectrum ”’ 
measured, the limb plates were marked perpendicular to the dispersion at 
one millimetre intervals from the inner edge of the spectrum. ‘The positions 
of each of the measured ‘‘spectra’’ were referred to these marks, and the 
position of the limb was determined visually with respect to them at both ends 
of the plate. ‘The image radius having been determined for the times of 
observation, a p/p, was calculated for each set of measurements. 

The measurements were made in the same manner as for the centre plates 
and mean positions of the seven lines were determined at each value of p/pp. 
Using the observed oxygen positions, these were corrected to the standard 
dispersion, assuming a linear correction factor. The difference between the 
corrected iron line positions and the standard iron line positions then gave the 
limb effect plus relative velocity corrections for the corresponding value of p/pp. 
The shifts so found for the three lines were then converted to equivalent Doppler 
velocities and averaged to give a mean shift for that value of p/py. 

I.5. Accuracy of measurement.—It was found that the velocity shifts determined 
as a function of p/p, from a single plate exhibited large fluctuations with amplitudes 
up to 0-2 km/sec. From an analysis of the micrometer settings in the centre 
spectra, it was concluded that the uncertainty in an individual reading of the 
position of an oxygen line was + 0-05 km/sec, and for an iron line + 0-07 km/sec. 
Similarly, an analysis of the limb measurements yielded an uncertainty in an 
iron line setting of +0-o6km/sec. From these figures, it is found that the 
uncertainty in the mean shift of all three solar lines at a given disk position due 
to measurement alone is + 0°05 km/sec. Since this uncertainty is much smaller 
than the observed fluctuations, the latter are attributed to the effects of local 
velocity fields on the Sun (15). 

From an analysis of micrometer measurements of the limb position with 
respect to the millimetre marks, it was concluded that the value of p/p) are 
uncertain by about +0-0003. At the limb, this corresponds to + 0-02 in cos @. 

1.6. Relative velocity corrections.—The corrections to be applied to the mean 
shifts obtained for each value of p/p, due to the motion of the observer and solar 
rotation were calculated using the Oxford Computing Laboratory ‘‘ Mercury”’ 


| 
| 
| 
| 
ote 
‘ 
’ 


No. 5, 1960 The solar red-shift 425 


computer. ‘lwo programmes were written for this purpose. The first provides 
the coordinates of the centre of the disk and various other solar parameters at 
the time of exposure. The input for this consists of the times of exposure and 
certain data directly as tabulated in the Nautical Almanac. The second 
programme uses the output of the first plus the observational p/p) values as input 
data and calculates for each observed position: cos@, solar latitude, correction 
for observer’s motion, correction factor for solar rotation. The calculations of the 
corrections follow the procedure given by H. H. Plaskett (18) and, following 
M. G. Adam (5), the law of solar rotation is taken to be 


R, = R, cos B(1 — 0-229 B) 


where R, is the solar equatorial velocity. In the calculation of the solar 
rotation corrections, R, was regarded as an undetermined parameter, to be 
supplied from the observations themselves. 

The mean shifts at each disk position were then corrected for the observer’s 
motion, giving at each point a residual shift 


A (8) = Ro. 


where R,.f(@) is the correction for solar rotation, f(#) being the calculated 
factor, and A(@) is the limb effect. This gives essentially 249 equations in 
R, and the A(#). The 1948 limb-effect curve showed, however, that A(6) 
varies only slowly with @ except very near the limb. Thus the observations were 
split into thirteen groups, each covering a small range in cos# over which A(6) 
could be considered approximately constant. For each group of observations, 
Ry and A(@) were then determined by a least squares solution. As in the 1959 
determination, the values of R, were fairly constant for the inner groups but 
decreased rapidly as the limb was approached. This was attributed to the 
effect of scattered light (5). From the four inner groups of observations, little 
affected by scattering, a mean value of 1-922km/sec was determined for Rp. 
Using this, the residual shifts were corrected for solar rotation, leaving values 
of A(8). 

1.7. Scattering corrections.—Since the positions of the observed points near 
the limb were uncertain by +0-02 in cos@, and scattering effects depend 
primarily on the limb being observed (east or west), the shifts obtained from the 
three east limb plates were averaged together in intervals of 0-02 in cos, and 
similarly for the west limb. These results are shown in the top half of Fig. 1. 
The large scatter in the observed mean shifts is attributed to local velocity-field 
effects (15), and it is seen that the west limb values are lower near the limb than 
those for the east limb. This is as would be expected, since scattering effects 
are far greater for the west limb where solar rotation and true limb effect 
combine to give a rapid variation of wave-length near the limb. 

The scattering effect was assumed to be represented mathematically by the 
simplified expression 


| "Hee dx 


+1 
[7 de 


u(x) = 


426 L. A. Higgs Vol. 121 
where u(x) is the true velocity shift (in the line of sight) at p/p) = |x|, (x) is the 
intensity at x, S(x—x,) is the scattering function and u(x) is the observed 
velocity shift. The limb darkening curve, /(x), used by M. G. Adam was chosen, 
and, in the absence of any other “‘ scattered light ’’ exposures of the sky spectrum, 
the scattering function used by her was adopted (5). 

To compute corrections for scattering, mean observed velocity-shift curves 
were obtained for the two limbs. Each was then smoothed and a true velocity 
shift curve was approximated by iteration, using the above expression for the 
scattering. From this, a scattering correction factor as a function of cos @ was 
determined for each limb. Using this, the two observed limb-etfect plots were 
then corrected, and the final values for the two limbs are given in the lower 
half of Fig. 1. Due to limitations in the iteration process, reliable scattering 
corrections could not be found for points nearer the limb than cos@=0-05, 
and hence the limb-effect results terminate at this point. 

1.8. Final limb effect.—The east and west limb results were then averaged 
together to give a final limb-effect curve. This is given in Fig. 2. It is evident 
that local velocity-field effects still persist in the final curve despite the averaging 
of six plates. A smooth curve (a parabolic formula in cos 6 fitted by least squares) 
through the observed points gives the values tabulated in Table III as a function 
of p/p). Included in this table are the corresponding values obtained by 
M. G. Adam for the solar equatorial region (§). ‘The new determination appears 


Final limb effect 
pipo A (km/sec) 
L. A. Higgs M. G. Adam 

0-998 0584 0°474 
0'995 0°524 
0'990 o'461 0°407 
O°417 
o'980 0°383 0°358 
0°975 0°355 0°339 
0°970 0°331 O°321 
0°263 
0'925 o'214 
o-18s 
o'875 o'168 
0-000 0°000 


to give larger shifts than those observed by her, but the difference between the 
two is entirely accounted for by differences in the standard solar positions used 
in the two determinations (see Table VI, Part III). The uncertainty in these 
standard positions precludes the measurement of any general difference 
<o°07 km/sec. Thus the latitude effect, if such exists, is of too small a 
magnitude to be detectable by these measurements. An analysis of the velocity 
shifts obtained from individual lines indicated small systematic differences 
between the shifts of the three lines, but these could also be accounted for by 
slight errors in the standard positions of the three lines. The question of the 
total red-shift at the limb, especially the possibility of a super-gravity shift, 


H 
= 
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cannot be answered until the red-shifts of these lines at the centre of the disk 
are known (Part IT). 
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Fic. 1.—Mean observed limb effect before and after corrections for scattering. 
Open circles: west limb. 
Filled circles: east limb. 
Continuous curve: final limb effect. 


Part II 


Determination of absolute wave-lengths and red-shifts 


II.1. Method.—The determination of red-shifts requires the absolute wave- 
lengths of both the solar lines and the corresponding terrestrial lines, formed in 
physical conditions as near as those existing in the solar photosphere as possible, 
e.g. vacuum arc. At Oxford, the method of ‘‘circular channels’’ has been 
developed as a powerful interferometric technique for wave-length determinations, 
the theory and applications of which have been presented in the literature 


(21, 2, 3, 4). 
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Essentially the method consists of using a spectrograph in conjunction with 
a Fabry-Perot interferometer, the spectrograph serving only as a monochromator 
for the interferometer. If the light being analysed has a continuous spectrum 
with absorption lines, the observed interferometric pattern consists of a set 
of circular ‘‘ heterochromatic’’ fringes, due to the continuous spectrum, which 
are crossed by absorption arcs forming part of another set of concentric circles 


Cosine 6 


Fic. 2.—Final limb effect. 
Open circles: mean of east and west limbs. 
Continuous curve: smooth limb effect. 


intersecting the first. In addition, the interference pattern of the cadmium 
red line (A 6438-4696) is superimposed on this pattern using a beam splitting 
cube. Overlapping is prevented by a system of diaphragms. For the measure- 
ment of arc wave-lengths, in which one is concerned only with a discrete line 
spectrum, the ‘‘heterochromatic’’ fringes do not exist and emission arcs now 
replace the absorption arcs. 

The interference patterns are measured on a specially adapted Hilger photo- 
measuring micrometer with a rotating plate stage. From detailed calculations, 
described in earlier papers, using the measured arc positions, the wave-lengths 
of the absorption or emission lines are compared with the wave-length of the 
cadmium line. Although complicated in its early stages, the method of circular 
channels has now become a standard procedure, and with the development of 
computer reduction of plate measurements is now a fairly rapid method of 
wave-length determination, 
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Il.2. Reduction of measurements.—The formerly tedious reduction procedure 
has been superseded by automatic reduction using the ‘‘ Mercury’’ electronic 
computer of the Oxford Computing Laboratory. A programme for this purpose 
has been written and is suitable for most conceivable wave-length determinations. 

The input data consist of general plate information and the plate measurements. 
The former include: the number of lines that have been measured; chromatic 
phase changes, camera focal length corrections and approximate wave-lengths 
for each line; and special instruction parameters to indicate the specific 
reduction procedure to be followed. The output consists among other things 
of a tabulation of the fractional orders of interference of all the absorption or 
emission lines and the absolute wave-lengths of all the lines. In addition, the 
mean error due to the scatter of the measurements is tabulated for each value. 
The total reduction time for a typical plate is reduced by this procedure to 
approximately twenty minutes, almost all of which is taken in preparing the 
input data on punched tape. 

11.3. Procedure.—The smaller Oxford telescope, giving an 18cm _ solar 
image, was used for the determination of absolute wave-lengths at the centre 
of the disk. The Sun was imaged on the spectrograph slit, the slit height being 
1°8 cm, and filters were used to equate the intensities of the solar and cadmium 
light beams passing through the interferometer. Two different spacers were 
used in the interferometer giving plate separations of 2-5 and 5mm. Fused 
silica plates were used and had been recoated at the beginning of the wave-length 
programme to give 7 per cent transmission and go per cent reflection in the red 
region. Seventy-four solar plates were taken, using Ilford panchromatic R40 
plates, at exposures varying from one to three minutes, depending on the spacer 
used. Of these, twenty-one plates taken on four different days were measured 
and reduced. Phase-change corrections and camera focal length corrections 
were adopted from earlier measurements (6). Both of these were small due to 
the small differences between the wave-lengths of the lines in question and the 
standard cadmium wave-length. For this reason, corrections for the dispersion 
of air in the interferometer were not necessary. 

For the arc wave-lengths, a vacuum arc run at about g cm Hg pressure replaced 
the coudé beam from the telescope. Mild steel pole-pieces were used and a 1omm 
arc was maintained at a current averaging 8 amperes. ‘The usual method of 
guiding the arc and focusing it on the slit was employed (2), only the central 
portion of the arc actually being used in order to eliminate any pole effects. Due 
to the low intensity of the arc lines, exposures of up to two hours were required 
and the cadmium intensity was correspondingly reduced using filters. Fifteen 
plates were taken, using Ilford HP3 plates, and twelve of these were measured 
and reduced. The corrections for phase change, etc., were the same as for the 
solar plates. The calculated absolute wave-lengths were reduced to absolute 
vacuum wave-lengths using the empirical pressure depression of spectral terms 
given by Babcock (g). The largest pressure correction for the three iron lines 
amounted to only o-g mA. 

11.4. Solar and arc wave-lengths.—The calculated solar wave-lengths 
obtained from the centre plates were corrected for the Doppler shifts due to the 
motion of the observer using the computer programme developed for Part I. 
Since some plates were of better quality than others, it was decided to weight 
each of the calculated wave-lengths on each plate according to the inverse square 
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of its calculated mean error for that plate. A similar method was used to determine 
the average arc-wave-lengths. ‘There were no serious differences in the wave- 
lengths determined considering the two spacers separately. For the solar 
wave-lengths, the 5 mm spacer gave wave-lengths o-6 mA longer than the 2-5 mm 
spacer, while for the arc wave-lengths, the 2-5 mm spacer values were 1-omA 
longer. The final absolute wave-lengths are given in Table IV, along with their 
standard deviations and the number of plates used in each determination. For 
comparison purposes, the I.A.U. wave-lengths are given for the solar lines and 
the Allegheny Observatory wave-lengths (10) of the arc lines are given. 

Il.5. Integrated disk wave-lengths.—Since it is known that local velocity fields 
may greatly affect wave-length determination on the Sun, the question arises 
whether the twenty-one solar plates were enough to eliminate these effects. 
As a check, wave-length plates were taken using integrated sunlight from all 
parts of the disk. Ideal conditions for integrated-light wave-length measurements 
would be attained if every part of the slit received truly integrated light and if 
every part of the spectrograph collimator received integrated light. The first 
condition was satisfied by removing the Cassegrain mirror of the telescope and 
allowing a parallel beam of light to fall directly on the spectrograph slit. The 
second condition, on the other hand, could not be satisfied. ‘To counteract this, 
an image rotator was placed before the slit and the circular channel plates were 
taken in ‘‘ mirror pairs’’. For the second plate of each pair, the image used for 
the first was inverted about the optical axis. In this way, wave-length shifts 
due to variations in the transmission properties of the grating across its surface 
were expected to be averaged out. Two pairs of plates were measured. 


TABLE IV 
Absolute wave-lengths and red-shifts 


Solar spectrum wave-lengths 


I.A.U, Observed No. Plates 
O, 6295°960 A 5°9603 + 00002 10 
Fe 6297°799 7°7989 + 0'0004 21 
O, 6299°228 9°2283 + 0'0002 21 
Fe 6301°508 15080 + 0°0003 21 
O, (6302-000) 19998 + 0°0001 21 
Fe 6302°499 2°5000 + 0°0003 21 


6302°764 2°7640 + 0°0001 10 


Iron arc wave-lengths 


Alleg. Obs. Observed No. .Plates Redeshift 
6297°790A 7°7918 + 0:0006 10 71t0°7mA 
6301°501 1°5008 + 0'0003 12 72+0°4 
6302°493 2°4938 + 0°0004 12 6-2+0°5 


The average wave-length values determined for each of these mirror pairs 
agreed well with each other, although there were differences (~ 4 mA) between 
the wave-lengths determined from the four plates individually. Since limb effect 
is known to concern only light coming from the extreme outer portion of the 
solar disk, one would expect the integrated disk wave-lengths to agree substantially 
with those determined for the centre of the disk. In fact, the agreement between 
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the mean wave-lengths from the four integrated disk plates (corrected for the 
observer’s motion) and the absolute wave-lengths determined at the centre of 
the disk is remarkable. Table V gives this comparison. This close agreement 
would seem to indicate that the absolute solar wave-lengths found at the centre 
of the disk are free from effects due to local velocity fields. 

11.6. Red-shifts—From the above absolute solar and arc wave-lengths, the 
red-shifts of the three lines may be determined. These are presented in the last 
column of Table 1V. These observed shifts are found to be slightly greater 
than those expected according to the empirical shift vs. line strength relation 
given by Nichols and Clube (16). 


TABLE V 
Comparison of solar wave-lengths 


Element Centre Integrated disk 
6295°9603 + 0:0002 A 59605 + 0'0002 
6297°7989 + 0°0004 7°7980 
6299°2283 +.0°0002 7°2288 + 0°0004 
6301°5080 + 0:0003 1°5086 + 
6301°9998 + 0°0001 2'0002 + 0°0003 
6302°5000 + 0:0003 2°4999 + 00007 
6302°7640 + 00001 2°7636 + 0°0003 


The results of Part I may now be extended to give the red-shift of the three 
lines across the disk (within the reservations to be indicated in Part IIT). 
A final limb-effect curve on an absolute scale is given in Fig. 3, in which the 
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Fic. 3.—Limb effect and the gravity shift. The uncertainty in the limb-effect determinations 
is indicated by the shaded band. 


Continuous curve: 1960 limb effect. 
Dashed curve: 1959 limb effect. 
Horizontal line: gravity shift. 
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1959 curve and the new determination are presented with their regions of 
experimental error. Also shown is the red-shift required by the principle of 
equivalence. From this diagram, it is evident that, allowing for experimental 
error, the shift observed at the extreme limb exceeds the gravity-shift, thus 
apparently ruling out the radial current interpretation of limb effect. 


Part Ill 


Line profiles 

III.1. Centre wave-lengths.—-In Part | of this investigation, standard micro- 
meter positions for a centre spectrum were adopted for the seven lines used, and 
in Part II, the determination of the absolute wave-lengths of these lines was 
presented. Using the absolute oxygen wave-lengths and micrometer positions, 
it was possible to calculate the micrometer positions corresponding to the absolute 
wave-lengths of the iron lines. It was found that the calculated micrometer 
positions of the three solar lines were, on the average, 0-005 mm to the red side 
of the adopted standard positions. ‘This corresponds to a wave-length difference 
of 4mA. The calculated and standard micrometer positions are compared in 
Table VI, in which the standard positions used by M. G. Adam (5) are also 
included. 


Taste VI 
Micrometer positions 


El. Standard positions Calculated positions 
M. G. Adam L. A. Higgs 
6296°0 A O, 0°000 mm mm mm 
6297°8 Fe 2°57! 2°570 275748 
6299°2 O; 4°576 4°576 4°5764 
6301°5 Fe 7°764 7°761 7°7683 
6302°2 O, 8-456 8°457 8-4569 
6302°5 Fe 9°154 97153 9°1572 
6302°8 O, 9°527 9°527 9°5268 


On account of the small portion of the spectrograph slit (corresponding to 
600 km on the surface of the Sun) used in the determination of the standard 
micrometer positions, one would expect some random fluctuations due to local 
velocity fields. Results obtained by A. B. Hart (15) indicate that r.m.s. velocity 
fluctuations of about 0-15 km/sec can be expected at a given point on the solar 
disk. Since ten plates were used in the determination of the standard positions, 
one might expect fluctuations in the mean positions of about +0001 mm. 
Since the discrepancy in question was five times this value, it seemed unlikely 
that it was due to unusually strong velocity fields at the centre of the disk when 
the ten spectra were taken. 

On the other hand, since the slit height used in the absolute wave-length 
measurements corresponded to 140000km on the surface of the Sun, one 
would expect all local velocity fluctuations to have been eliminated. The close 
agreement between the wave-lengths determined at the centre of the disk and 
those obtained using integrated light from the whole disk tended to confirm this. 


It thus seemed difficult to accept the possibility that the absolute wave-lengths 
were in error. 
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Since an asymmetric apparatus function of a spectograph can produce 
wave-length shifts of spectral lines (17), the possibility of such an effect in the 
grating spectra was considered. Using the Oxford Moll-type microphotometer, 
sixteen tracings of typical spectra were taken, in the forward and reverse directions, 
both at the centre of the disk and at several points near the limb. The dimensions 
of the microphotometer slit were chosen so as to correspond to the dimensions 
of the eyepiece diaphragm used in the micrometer measures. The tracings of 
a typical centre spectrum showed that the oxygen lines were approximately 
symmetric while the solar lines were definitely asymmetric, having pronounced 
violet flanks. Since this asymmetry is of such a magnitude as to be scarcely 
detectable in a diagram of the actual line profile, the following method of 
presentation will be used. The profile of a line is represented by a curve giving 
the mid-point of the profile at each value of the incident intensity, J. This is 
illustrated by the curve AB in Fig. 4(a). In this manner, the profile at the 
centre of the strongest iron line is presented in Fig. 4 (6). 

To determine an approximate apparatus function for the 30000 line/inch 
Bausch and Lomb grating, it was decided to use the terrestrial oxygen lines as 
essentially monochromatic absorptions. The apparatus function is then given 
by the observed oxygen profiles. Allen (8), who has used atmospheric oxygen 
lines for similar purposes, gives theoretical widths for these lines of about 18 mA, 
which indicates the degree of approximation involved. Only the limb tracings 
were used for this purpose because the centre spectra were too over-exposed 
to show any possible asymmetry in the oxygen lines. All profiles were reduced to 
true intensity profiles using a characteristic curve for R50 plates, obtained in 
another investigation. From profiles of the oxygen lines A6296-0 and A6299°2 an 
approximate apparatus function was obtained. Having a half-width of 69mA, 
it is symmetric to a width of about 150 mA, after which it develops a violet wing. 
The uncertainty in the determination of the wings is about 1 per cent of the 
central intensity, whereas the asymmetry in the wings is three times this figure. 
‘Thus a true asymmetry is indicated. 

Using this apparatus function, it was shown that the observed asymmetry 
of the solar lines at the centre could be accounted for by purely instrumental 
effects, the true solar profiles being symmetric. In addition, it was estimated that 
the above apparatus function would shift the part of the solar lines used in 
micrometer settings 3-4mA to the violet from their true wave-lengths. This 
explains the apparent discrepancy between the grating and interferometric 
wave-lengths. Since this instrumental shift of grating wave-lengths is 
independent of the point of observation on the solar disk, it has no effect on 
micrometer limb-effect determinations. 

III.2. Limb profiles.—To investigate the profiles of the solar lines near the 
limb, tracings of six limb spectra were measured, each at different values of 
p/py. ‘The east-limb plates were used for these in order to reduce the effects of 
scattered light. Fig. 4(c) to (f) shows the observed “‘ profile’’ of the strongest 
iron line (A6301°5) as a function of p/p). There seems to be a systematic change 
in profile as one approaches the limb, with the asymmetry reversing in the core 
of the line, indicating effects other than instrumental slit pattern. Although 
these curves all refer to the strongest iron line, the other two iron lines show 
similar behaviour. 
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Several causes of the difference between centre and limb profiles were 
considered. The most plausible of these was that the scattering of light might 
create very asymmetric line profiles at the limb where both the intensities and 
wave-lengths of the lines are varying rapidly. To test this, an ‘‘ observed’’ 
profile was computed for p/p,=0-997, assuming the true line profile to be 
symmetric. In this calculation, the same scattering function was used as in 
Part I, the apparatus function was taken to be that above, and the wave-lengths 
were calculated using the limb-effect curve in Part I. It was found that scattering 
could not account for the observed asymmetry of the lines at the limb, its only 
effect being a slight shift of the whole line to the violet. In a similar manner, it 
was shown that the finite height of the microphotometer slit could not account 
for the observed profiles. 
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Fic. 4.—Measured asymmetry of the 6301°5 A iron line. 
(a) Explanatory profile. 
(6) Centre. 
(c) p/pp=0'988. 
(a) =0°994. 


It thus appeared possible that the solar lines at the extreme limb are formed 
with asymmetric profiles. The unexplained asymmetry is seen to set in at a 
value of p/p, of about 0-98 (about 20” from the limb). Such an asymmetry has 
never been previously observed* ; on the other hand, observations so near the 
limb are few. Assuming that such true asymmetry existed, an attempt was made 
to correct the observed profiles for the apparatus function and thus determine 
the amount of the true asymmetry. Unfortunately, the asymmetries concerned 
were too small and the apparatus function insufficiently well determined for any 


* Although asymmetry of infra-red oxygen lines has been detected at the centre by H. H. Voigt, 
Zs. f. Ap., 4, 365, 1956. 
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iteration procedure to yield reliable results. It did appear, however, as if the 
solar lines at the limb have very pronounced red flanks. 

The asymmetry observed in limb lines, if real, completely changes our views 
as to the nature of limb effect. A wave-length shift to the red will in all probability 
be produced by the strong red flanks which are indicated for the limb lines, but 
the interest is now changed from wave-length measurement, an unsuitable 
means of measuring line asymmetry, to a photometric survey of the nature and 
origin of the asymmetry observed. To this end, observations are at present in 
progress with a Babcock grating, the apparatus function of which is being 
accurately determined, giving greater dispersion and resolving power than the 
Bausch and Lomb grating so far employed. 
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Summary 


A survey has been made of the orbits of faint sporadic meteors. ‘lhe 
results show that, in addition to meteors moving in orbits similar to those of 
the comets, there are many in orbits of very short period. ‘The more eccentric 
of these orbits are concentrated towards the plane of the ecliptic, but some 
are nearly circular and have inclinations near 60°. "The differences between 
these results and those previously published for bright meteors are dis- 
cussed, and a qualitative explanation is given in terms of the action on near- 
parabolic orbits of solar radiation and the major planets. 


1. Introduction 

The orbit of a meteor may be calculated if both its velocity and direction of 
motion are measured as it enters the Earth’s atmosphere. Measurements of 
these quantities have been made visually, but it is now generally recognised that 
visual estimates of meteor velocity are not accurate enough to be used to derive 
the orbit. Except for a single orbit measured by McKinley and Millman (1) 
by a radio-echo method applicable to only a few of the brightest meteors, all the 
reliable measurements of the orbits of individual meteors so far published have 
been made using multiple station photography (see, for example, Whipple (2)). 
The photographic method has great accuracy, but its use in the study of faint 
meteors is difficult ; thus the §1 sporadic meteors for which Whipple has published 
orbits were all brighter than visual magnitude 0. Although cameras able to 
record meteors of magnitude +4 have now been developed (3), there is no 
immediate prospect of photographing meteors in the magnitude range +6 to 
+9. Radio-echo observations of meteors as faint as this have been made but, 
because no means of measuring the direction of motion of an individual meteor 
has been available, they have previously been limited to statistical surveys of 
meteor velocities (e.g. Almond, Davies and Lovell (4), McKinley (§)) and of 
radiants (Hawkins (6)). 

We have developed a method of using a transmitter and three spaced receivers 
to measure both the velocity and direction of motion of a meteor. This method 
has now been used in a survey of the orbits of faint meteors (visual magnitude 
approximately +7), the results of which are presented below. 

2. The method and apparatus 

As the apparatus used to make the observations, and the method of 
calculating a meteor orbit from them, have been described in a previous paper 
(7), only a brief account will be given here. 

The apparatus consists of a pulsed transmitter, together with three receivers 
which are situated on the ground approximately at the vertices of an isosceles 
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right-angledtriangle. Ifa meteor trail is formed in such a position that the radiation 
from the transmitter is reflected specularly from it to the receivers, the received 
echoes may show Fresnel diffraction patterns. At any two stations these patterns 
are separated in time by an interval equal to that taken by the meteor to travel 
between the points of specular reflection. The product of this time interval 
with the meteor velocity, which itself may be found from the range of the echo 
and the time-scale of the diffraction patterns, gives the distance along the trail 
between the points of reflection to the two stations. The ratio of this distance 
to that separating the receivers on the ground is equal to half the cosine of the 
angle between the directions of the trail and of the line joining the stations. 
From the three diffraction patterns are found two independent direction cosines, 
defining the direction of motion of the meteor. As the velocity of the meteor and 
the time of its occurrence are known, the orbit may then be computed. 

The transmitter used in the survey radiated at 36:3 Mc/s, emitting 600 
pulses a second with a peak power of 1tookw. The transmitting aerial had 
an absolute power gain of 28, a beam width of + 12°, and was directed slightly 
north of east at an elevation of 30°. An aerial of gain 12 in this direction was 
used at each of the three receiving stations. One receiving station was near 
the transmitter; the other two were approximately south and east of this at 
distances of about 3°5 kilometres. Auxiliary transmitters relayed to the first 
station the signals received at the others, and echo pulses received at all 
three stations were displayed on adjacent cathode-ray tubes and recorded 
photographically. 

The orbits measured in the survey were calculated, following the method 
of Porter (8), using the Manchester University digital computer. 

3. Experimental details 

(a) Times of observation.—The survey was made during a twelve-month 
period, using the apparatus for 24 hours each month. Details are given in 
Table I. 

| 
Year and month Observation period (U.T.) No, of orbits computed 
‘Time of start Time of finish 
days hours mins. hours mins. 
1954 May 20 00 20 247 
June 22 
2210 22 
July 27. «93 28 305* 
August 26 10 27 206 
September 26 27 234 
October 21 i 22 278 
November 23 
23 23 166 
24 ol 24 
December 20 «10 21 159 
1955 January 24 10 25 128 
February 20 21 140 
March 21 12 22 1607 
April a2 (oF 22 168 
* Plus 35 5-Aquarid orbits. 
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(6) The visual magnitudes and masses of the meteors.—A knowledge of the 
visual magnitudes of the meteors observed in the survey is useful when comparing 
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the results with those obtained photographically. The magnitudes may be 
estimated from the electron densities in the trails, by comparison with previously 
published work on the relation between luminosity and ionization in meteor trails. 

(i) The electron densities in the trails.—The amplitude and range of the radio 
echo received from a meteor trail are related to the line density of electrons in 
it by the formula of Lovell and Clegg (9); 


a=2V3nA 


Here « is the number of electrons per cm of trail at the reflecting point, 
A is the echo amplitude (in terms of the r.m.s. noise at the receiver input), 
R is the echo range, k is Boltzmann’s constant, T the receiver equivalent input 
temperature, B the bandwidth, P the power radiated during the transmitter 
pulse, A the wavelength, and G,, G, are respectively the absolute power gains 
of the transmitting and receiving aerials in the direction of the reflecting point 
on the trail. 

Of these quantities, P, A and B (150ke/s) are known parameters of the 
apparatus. 

The equivalent temperature T was about 6000°K. Most of the noise was 
extra-terrestrial in origin, and T varied as different parts of the galaxy were 
viewed by the receiving aerials. With the wide-beam aerials used, however, 
the variation was not large and has been ignored. 

The gains G,, G, were not known accurately because, although the charac- 
teristics of the aerials are calculable, the directions of the echoing points were 
not measured. The elevation » was however estimated from the range R of 
the echo, since siny ~h/R and the height of the echoing point could be assumed 
to be about g95km. As the gain is unlikely to have been much less than the 
maximum possible for elevation , a fairly good estimate of G, G, could be 
made. 

An approximate value of the electron density in the trail at the reflecting 
point was thus obtained. For most of the meteors observed in the survey it was 
between 5 x 10'°and 2 x 10" percm. The maximum trail density, «max, usually 
exceeded that measured, and it can be shown from the theory of Herlofson 
(10) that the average value of &max/« is likely to be about 1-4; from this it has 
been estimated that, for about three-quarters of the meteors observed, «max 
was between 5 x 10" and 4x 10" per cm. 

(ii) The relation between visual magnitude and electron density.—The 
relationship between the visual magnitude of a meteor and the number of 
electrons in its trail has been studied by Kaiser (11). Using the results 
of Millman, McKinley and Burland (12) for Perseid meteors, he concluded 
that a 5th or 6th magnitude meteor would produce a maximum of 10" electrons 
per cm of trail. Densities of 5 x 10’ and 4 x 10" per cm therefore correspond to 
Perscids of visual magnitudes of +8 and +6 respectively. The Perseids 
have velocities near 60km/sec. According to Kaiser, the brightness of a 
meteor may be expected to be proportional to kv*8~xmax where v is its velocity, 
k the fraction of its kinetic energy converted into light, and f the probability 
that an evaporated meteor atom will give rise to a free electron, ‘The extra- 
polation of the Perseid results to meteors of other velocities thus depends on 
the variation with velocity of kv*8~'. Kaiser finds that, if 8 varies as v”, then 
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right-angledtriangle. Ifa meteor trail is formed in such a position that the radiation 
from the transmitter is reflected specularly from it to the receivers, the received 
echoes may show Fresnel diffraction patterns. At any two stations these patterns 
are separated in time by aninterval equal to that taken by the meteor to travel 
between the points of specular reflection. The product of this time interval 
with the meteor velocity, which itself may be found from the range of the echo 
and the time-scale of the diffraction patterns, gives the distance along the trail 
between the points of reflection to the two stations. The ratio of this distance 
to that separating the receivers on the ground is equal to half the cosine of the 
angle between the directions of the trail and of the line joining the stations. 
From the three diffraction patterns are found two independent direction cosines, 
defining the direction of motion of the meteor. As the velocity of the meteor and 
the time of its occurrence are known, the orbit may then be computed. 

The transmitter used in the survey radiated at 36°3 Mc/s, emitting 600 
pulses a second with a peak power of 1tookw. The transmitting aerial had 
an absolute power gain of 28, a beam width of + 12°, and was directed slightly 
north of east at an elevation of 30°. An aerial of gain 12 in this direction was 
used at each of the three receiving stations. One receiving station was near 
the transmitter; the other two were approximately south and east of this at 
distances of about 3°5 kilometres. Auxiliary transmitters relayed to the first 
station the signals received at the others, and echo pulses received at all 
three stations were displayed on adjacent cathode-ray tubes and recorded 
photographically. 

The orbits measured in the survey were calculated, following the method 
of Porter (8), using the Manchester University digital computer. 


3. Experimental details 

(a) Times of observation—The survey was made during a twelve-month 
period, using the apparatus for 24 hours each month. Details are given in 
Table I. 


I 
Year and month _ Observation period (U.'T.) No, of orbits computed 
‘Time of start ‘Time of finish 
days hours mins. days hours mins. 
1954 May 20 00 20 247 
June a8 23 22 
aa. 22 
July 27 «(03 28 
August 27 206 
October 21 il 22 278 
November 23 
23 (04 23 166 
24 O1 24 
December 20. 21 159 
1955 January 24 25 128 
February 20 21 140 
March 21 2 22 167 
April 22 7 22 168 
* Plus 35 5-Aquarid orbits. 


(6) The visual magnitudes and masses of the meteors.—A knowledge of the 
visual magnitudes of the meteors observed in the survey is useful when comparing 
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the results with those obtained photographically. The magnitudes may be 
estimated from the electron densities in the trails, by comparison with previously 
published work on the relation between luminosity and ionization in meteor trails. 

(i) The electron densities in the trails. —The amplitude and range of the radio 
echo received from a meteor trail are related to the line density of electrons in 
it by the formula of Lovell and Clegg (9); 


RRTB 
x= 2V30A( ) 


Here « is the number of electrons per cm of trail at the reflecting point, 
A is the echo amplitude (in terms of the r.m.s. noise at the receiver input), 
R is the echo range, k is Boltzmann’s constant, T the receiver equivalent input 
temperature, B the bandwidth, P the power radiated during the transmitter 
pulse, A the wavelength, and G,, G, are respectively the absolute power gains 
of the transmitting and receiving aerials in the direction of the reflecting point 
on the trail. 

Of these quantities, P, A and B (150kc/s) are known parameters of the 
apparatus. 

The equivalent temperature 7 was about 6000°K. Most of the noise was 
extra-terrestrial in origin, and T varied as different parts of the galaxy were 
viewed by the receiving aerials. With the wide-beam aerials used, however, 
the variation was not large and has been ignored. . 

The gains G,, G, were not known accurately because, although the charac- 
teristics of the aerials are calculable, the directions of the echoing points were 
not measured. The elevation 7 was however estimated from the range R of 
the echo, since siny ~h/R and the height of the echoing point could be assumed 
to be about g5km. As the gain is unlikely to have been much less than the 
maximum possible for elevation 7, a fairly good estimate of G, G, could be 
made. 

An approximate value of the electron density in the trail at the reflecting 
point was thus obtained. For most of the meteors observed in the survey it was 
between 5 x 10'°and 2x 10" percm. ‘The maximum trail density, «max, usually 
exceeded that measured, and it can be shown from the theory of Herlofson 
(10) that the average value of amax/« is likely to be about 1-4; from this it has 
been estimated that, for about three-quarters of the meteors observed, «max 
was between 5 x 10” and 4x 10" per cm. 

(ii) The relation between visual magnitude and electron density.—The 
relationship between the visual magnitude of a meteor and the number of 
electrons in its trail has been studied by Kaiser (11). Using the results 
of Millman, McKinley and Burland (12) for Perseid meteors, he concluded 
that a 5th or 6th magnitude meteor would produce a maximum of 10” electrons 
per cm of trail. Densities of 5 x 10" and 4 x 10"! per cm therefore correspond to 
Perscids of visual magnitudes of +8 and +6 respectively. The Perseids 
have velocities near 60km/sec. According to Kaiser, the brightness of a 
meteor may be expected to be proportional to kv*8~'amax where v is its velocity, 
k the fraction of its kinetic energy converted into light, and 8 the probability 
that an evaporated meteor atom will give rise to a free electron. The extra- 
polation of the Perseid results to meteors of other velocities thus depends on 
the variation with velocity of kv°8~'. Kaiser finds that, if 8 varies as v", then 
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n=o+o6, If this result is combined with Jacchia’s assumption that & is 
proportional to , the brightness of a meteor trail of given electron density is 
found to be proportional to v*, 

More recently Whipple (13) has suggested that & varies as v*. From the 
relative abundances of radio and photographic meteors of different velocities 
he finds Bw", and therefore ke®8-' to be proportional to v. It is shown in 
this paper, however, that the published distributions of the velocities of radio 
meteors are affected by selection favouring the measurement of the faster 
meteors, and also that the orbital distributions of faint (radio) and bright 
(photographic) meteors are different. If these factors are allowed for, kv®f ' 
is found to depend on a higher power of v, 

‘Thus, although the evidence is hardly conclusive, it seems probable that 
ko" 8B! varies as a fairly high power of v. If the fourth power relation is 
assumed to be the correct one, a maximum electron density of ro" per cm, 
corresponding perhaps to a 7th magnitude Perseid meteor, could also result 
from a 2okm/sec meteor of visual magnitude +11°5. Some of the slower 
meteors observed in the survey may, therefore, be extremely faint visually. 

(iii) The masses and dimensions of the meteors.-The masses and radii of the 
meteors are also of interest, since solar radiation affects cach orbit to an extent 
which depends on the ratio of the surface area of the meteor to its mass. 

Kaiser (14), using the results of Millman and Jacchia, estimates a meteor 
producing 10" electrons per cm of trail to have a mass between 10? and 10 “gms. 
If, as he suggests, the probability of ionization is constant, this applies to meteors 
of any velocity. ‘The survey meteors should then have masses between 1o~4 
and 10° °gmsa, radii of the order of 10°-% ems, and a ratio of surface area to mass 
of perhaps socm® per gm. 

More recently, the Harvard group have shown from their photographic 
observations that the fainter metcors are probably loose agglomerations of 
particles and not, as previously supposed, dense spherical bodies of stone or 
iron, ‘They estimate average densities of the order of o-2zgm/cem", and 
masses perhaps tooo times as great as those previously accepted. The 
surface-area/mass ratio, however, is found directly from the luminous efficiency, 
and the previously accepted values may well be correct. 

(c) Instrumental and astronomical selection.—In order to tind the orbital 
distribution of the meteors intercepted by the Earth, the orbit of each meteor 
measured in the survey must be given a weight inversely proportional to the 
probability that the meteor, having entered the atmosphere, will produce a 
trail suitable for measurement, If the true distribution of meteors in the solar 
system is required, it is necessary to include a further weighting factor, inversely 
proportional to the probability that the meteor will collide with the Earth, 

‘These factors are referred to below as the instrumental and the astronomical 
weighting factors. 

(i) The instrumental weighting factor,-The orbit of a meteor can be cal- 
culated only when an echo showing a measurable Fresnel ditfraction pattern is 
recorded from each of the three stations. Echoes whose amplitudes exceed 
five times the receiver noise, and which show at least three or four zones of the 
diffraction pattern, are normally required. 

‘The amplitude of the echo depends on the position and direction of the 
trail, and on the electron density within it. Echoes are received only from 
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trails formed in the region illuminated by the transmitting aerial: i.e. the 
meteor must form a trail within the ‘collecting volume’ of the apparatus. 
The trail must be oriented during its formation so as to reflect radiation specularly 
to the receiving stations. This condition requires nearly all the meteors 
observed in the survey to come from a strip of the celestial sphere determined 
by the directional properties of the aerials. This strip stretched from the 
north to south, passing about 30° west of the zenith, and was about 20” wide. 
The electron density in the trail has to be sufficient to give an echo large enough 
to be measured, but not so great as to saturate the receiver. Because of this, 
most of the survey meteors produced trails having maximum electron density 
between 5 = 10" and 4» 10" perem, 

Very few orbits were measured from trails having more than 2:10!* electrons 
per cm. Such trails give rise to ‘long-duration’ echoes, which rarely show 
measurable diffraction patterns; the orbits obtained from them have not been 
included with those in the main survey, on account of the different magnitude 
range and selection factors involved. 

The echo amplitude, for a meteor producing a trail of given density, thus 
depends on the position of the radiant. A dependence on velocity may also 
arise. The fastest meteors, because they evaporate at the greatest heights, may, 
on account of the rapid diffusion of their trails, give echoes of reduced amplitude 
and in extreme cases it may be impossible to observe any diffraction pattern, 
However, for the survey meteors this dependence of echo amplitude on velocity 
is expected to be slight and is neglected. 

About three-quarters of those ‘ short-duration ' echoes which had amplitudes 
suthcient for measurement were nevertheless unsuitable for orbit determination 
because the diffraction patterns were cither distorted or badly defined. An 
investigation has been made of the causes of these imperfections in the observed 
patterns. It is believed that echoes showing distorted patterns come from 
trails formed in regions of atmospheric turbulence where the wind velocity 
changes non-uniformly with distance along the trail. Such echoes usually 
give inaccurate measures of meteor velocity and radiant, and were therefore 
rejected. Since the time taken for a meteor to produce the part of the trail 
which gives the diffraction pattern is inversely proportional to its velocity, the 
distortion, and therefore the chance of rejection, should be less for a fast meteor 
than for a slow one. As upper-atmosphere winds are known to be stratified 
horizontally, the distortion might also decrease with increasing distance of the 
meteor radiant from the zenith. Both these tendencies may be partially offset 
by the fact that meteors moving with high velocity, or with radiants near the 
horizon, produce trails at greater heights, in regions where turbulence tends to 
be greater. 

Echoes which show badly defined diffraction patterns are thought to arise 
from trails produced by fragmenting meteors. If a meteor disintegrates in 
the atmosphere into ionizing fragments, the smaller pieces are retarded more 
than the larger ones and a train of particles is formed. If the length of this 
train becomes comparable with, or exceeds, those of the early Fresnel zones of 
the trail, a diffraction pattern is unlikely to be seen, It can be shown that, if 
the individual fragments of the meteor themselves obey Herlofson’s equations 
(10), the length of the train of particles is likely to be proportional to (v* cos x)~"', 
where © is the meteor velocity and y is the zenith distance of its radiant. The 
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proportion of echoes showing no well-defined diffraction pattern should therefore 
be strongly dependent on velocity, and slightly dependent on radiant position. 
As in the case of the distorted patterns, the proportion of meteors rejected on 
account of fragmentation will be greatest for slow meteors. The meteors coming 
from near the zenith will in this case, however, have the best chance of measurement. 

Thus the quality of the diffraction pattern depends mainly on the velocity 
of the meteor, whereas the amplitude of the echo is a function chiefly of the 
radiant position. The probability that the orbit will be measured may therefore 
be written 

P(HA,5,v)xp,(HA, ).p,(v) 


where the meteor has velocity v and comes from a radiant at declination 6 
and hour angle HA, and p,, p, are functions of the radiant position and meteor 
velocity respectively. 

The survey was made by observing for several periods, each of one sidereal 
day. Since, during such a period, all radiants at a given declination have equal 
probabilities of observation, it is permissible to write, for sporadic meteors, 


P(5,v)<p,(5). 
The instrumental weighting factor, W(5,v) which is the inverse of P, may 
therefore be regarded as the product of two independent factors w(5) and w(v). 


The determination of w,: 

The factor is inversely proportional to the probability p, that a meteor of 
given velocity, coming from a radiant at declination 5 will be observed during 
a 24-hour period. It cannot be calculated reliably from the directivity of the 
aerial systems, because any dependence of diffraction pattern quality on radiant 
position is not then taken into account. The factor has been estimated, however, 
by making the following assumptions : 

(i) The sensitivity of the apparatus, as a function of radiant position, follows 
contours whose shape is independent of the meteor velocity. This is a sufficient 
condition for W to be regarded as the product of w, and w,, 

(ii) The contours are symmetrical about the vertical plane of symmetry of 
the aerial systems. ‘Thus for any area south of this plane (which lies in azimuths 
80° and 260°) there is an equal area of the same sensitivity in the northern 
half of the sky. 

(iii) The true distribution of sporadic radiants in any region of the celestial 
sphere does not vary as the region passes through the area to which the apparatus 
is sensitive. For any fixed declination, therefore, the observed density of radiants 
taken with respect to hour angle, is proportional to W and to the true density. 

Assumption (iii) is applicable only to radiants which have been corrected 
for the effects of the Earth’s attraction and rotation, whereas (i) does not neces- 
sarily hold when these radiants are considered, since the corrections are them- 
selves functions of meteor velocity. However, the corrections are usually 
small, except for meteors from radiants near the horizon, which are in any case 
rejected because the errors of radiant measurement are large. 

Assumptions (i), (ii) and (iii) may therefore be taken to apply to the 
distributions of corrected radiants. The factor w, has been found as follows: 

The corrected radiants of all the meteors for which orbits were found 
(except those measured during the incomplete run of 1954 June) were sorted 
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into 10° intervals of declination and hour angle. For each 10° strip in declina- 
tion, curves were drawn showing the distribution of radiants in hour angle, 


These are shown in Fig. 1. 


50 


— 


120° 60° 
Fic. 1.— Distributions in hour angle of observed radiants in different declination groups. 
Abscissa: hour angle in degrees. 

Ordinates: meteor numbers per 100 square degrees. 


The curve for the strip between declinations 40° and 50° (AB in Fig. 2) 
is derived from a large number of meteors and is a convenient starting-point 
for the calculation of w,. We define the sensitivity S(5,HA) of the apparatus 
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as being the number of meteors observed in unit time from a unit area in this 
strip. The unit of time is taken to be the total observing time in the survey, 
and the unit of area is the area of a 10° x 10° section of the strip. Curve (4) 
of Fig. 1 then represents the variation of S as a function of HA. 


Fic. 2.—Diagram of northern celestial hemisphere with curves of constant elevation. ZAD 
represents the plane of symmetry of the serial system, AB and AC are symmetrically placed ahout 
this plane. 


From this curve, and assuming (i) and (ii) above, the value of S is known also 
for the strip AC of Fig. 2, this being symmetrical about azimuth 260° with AB. 
Since at any fixed declination the number of meteors observed per unit area 
is proportional to S, the value of S may be found for all points in the 10° strips 
which are intersected by AC. By repeating this process, starting with strips 
other than AB, the value of S was found for every point from which meteors 
were recorded. 

The probability p, is then given by 


) 


the integral being taken over one complete cycle of hour angle. It was 
found for each declination by adding the calculated values of S for each 10° of 
hour-angle, and from its reciprocal the values of w, given in Table II were 
obtained. 
Il 
Declination Weight Declination Weight 
(degrees) (degrees) 
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— 20 —16 12°5 
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—12 —8 


HS 


Desh 

\ 

\ 
~ 
; 

MH 180 HA- 90 

4 

oe © to 4 5°6 52 to 56 

e 
4 tos 4°4 56 to 60 

8 to 12 3°5 60 to 64 
eae 12 to 16 2°7 64 to 68 

ey 16 to 20 2°3 68 to 72 

> 72 


No. 5, 1960 Measurements of the orbits of faint sporadic meteors 445 


The determination of w,: 

The weighting factor w, is inversely proportional to the probability p that 
a meteor of velocity v will give diffraction patterns suitable for orbit measure- 
ment. 

The factor was found from the records of 1954 May 20-21, by measuring 
the velocities of all the meteors which gave echoes of amplitude exceeding 
five times the receiver noise level. The velocities were measured from the best 
of the recorded diffraction patterns, the accuracy being perhaps 10 per cent. 
It was found that 94 per cent of the meteors provided at least one diffraction 
pattern giving the velocity to this degree of accuracy, and for a further 5 per 
cent a very rough estimate was possible from the initial rate of rise of echo 
amplitude. Only 1 per cent of the meteors gave no indication of their velocity. 

The percentage of meteors of various velocities for which orbits were 
measured is shown in Fig. 3. 


ao Su bo 7° 


Fic. 3.—The proportion of meteors accepted for orbital determination, as a function of velocity. 
Abscissa: observed velocity, km|sec. 
Ordinate: percentage of meteors used. 

The probability of orbit measurement was found to be approximately 
proportional to the velocity of the meteor. The weighting factor w, has 
therefore. been chosen to be inversely proportional to v. In calculating the 
various distributions given in Sections 5 and 6, a weight W=w,w, has been 
ascribed to each meteor. In the distributions of orbital parameters, which 
would be expected to show symmetry about the plane of the ecliptic, it has 
been found convenient to consider only those meteors whose geocentric radiants 
were north of the ecliptic. The data for the other hemisphere are inevitably 
incomplete, as meteors having radiants south of declination 22°S were not 
included in the survey, and the weighting factor w, is rather uncertain for 
meteors having radiants much south of the equator. 

The astronomical selection.—A weighting factor C, inversely proportional 
to the probability of collision of a particle with the Earth, has been given by 


Whipple (2): i 
I 


where P is the probability of collision per revolution of the particle, V,, is the 
geocentric velocity of the particle (neglecting the effect of the Earth’s attraction), 
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V,. is the relative velocity of the particle at the top of the Earth’s atmosphere, 
R is the radius of the Earth, and p=q(1 +e), where g is the perihelion distance 
(in astronomical units), e the eccentricity, a the semi-major axis and i the 
inclination of the particle’s orbit. 

This weighting factor may be applied to the meteor orbits intercepted by 
the Earth to obtain the true distribution of particles amongst these orbits. 


4. The experimental errors 

(i) Errors in the apparent radiants and velocities.—The accuracy with which 
the apparent radiant and velocity are measured with the survey apparatus has 
been discussed in a previous paper (7). It was shown that the errors are 
mainly due to atmospheric, and not instrumental causes. 

The errors in radiant position arise because the measured time-displace- 
ments between the diffraction patterns recorded at the three stations are affected 
when the trails are distorted by non-uniform winds. From observations of the 
Geminid shower of 1954 December, it was found that the standard deviation in 
measurements of the time-displacements is about +5%. The error in the 
radiant depends on position, but is typically +2°. 

The error in geocentric velocity arises from the uncertainty in the correction 
that is added to the measured velocity to allow for the retarding effect of the 
atmosphere. ‘This uncertainty was found to be + 2km/sec for Geminid meteors 
of velocity 36km/sec. Similar results have since been obtained from observa- 
tions on 1954 July 26—28 of the 5-Aquarid meteors, whose velocity is 43 km/sec. 

Although the errors have not yet been measured for meteors whose velocity 
differs greatly from 35 km/sec, it is thought that they decrease with increasing 
meteor velocity. We assume them to vary inversely as the meteor velocity, 
and allow for errors of 5 x 35/v per cent in the time displacements, and of 
2 x 35/v km/sec in the velocity before entering the atmosphere, for a meteor 
whose velocity is «km/sec. 

(ii) Errors in the calculated orbits.—-It would be dithcult to derive and use 
a general expression for the effect on the calculated orbits of given errors in the 
observed quantities. ‘The errors in the orbits may best be found separately 
for each meteor, by modifying appropriately the values of the measured velocity 
and time displacements and then re-calculating the orbit. ‘The results of some 
calculations of this type are shown in Table IIl. The standard deviations in 
the orbital parameters are given for a selection of meteors; the errors in the 
measured quantities are taken to be of the size assumed in (i) above. 

It is apparent from the table that the calculated values of 1/a and g may 
usually be relied upon to +o-1. As would be expected, the errors in g’ become 
large when near-parabolic orbits are considered. ‘The eccentricity e is usually 
measured to + 0°05, and the inclination to perhaps +10°. The argument of 
perihelion, w, can often be found to +10°; though for nearly circular orbits 
the error is likely to be much greater. However, with such orbits the ditference 
between the aphelion and perihelion distances is small, and w ceases to have 
much significance. 


5. The distribution of geocentric radiants and velocities 

Previous observations of faint sporadic meteors have been restricted to 
separate statistical surveys of radiants and of velocities. It is therefore of 
interest, before discussing the measured orbits, to consider the radiants and 
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velocities measured in the present survey, and to compare them with the earlier 
data. 

(a) The geocentric radiants. (i) The measured distribution.—The distribution of 
geocentric radiants is shown in Figs. 4, 5 and 6. The radiants have been 
corrected for the effects of the Earth’s attraction and rotation. So that the 
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Fic. 4a.—Observed distribution of radiants in ecliptic longitude for each month. 

Abscissa: longitude, measured from the apex of the Earth’s way. 

Ordinate: observed mumbers per 10° of longitude. 
radiants of meteors observed on different days could be included in one dis- 
tribution, the ecliptic latitude § of the radiant, and the differences in longitude 
A—A between the radiant and the apex of the Earth’s motion, have been chosen 
as co-ordinates. 

Fig. 4 shows separately the distributions obtained in the twelve periods of 

observation. In Fig. 5 the results of the different periods have been combined, 
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and the numbers of meteors having radiants at various latitudes are indicated. 
The complete latitude distribution is given in Fig. 6. In each figure the 
distribution is shown both before and after the instrumental weighting factors 
W, and W, have been included. 
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Fic. 46.—Distribution of radiants in ecliptic longitude for each month, corrected for observational 
selection. 
Abscissa: longitude measured from the apex of the Earth's way. 
Ordinate: percentage of total distribution per 10° of longitude. 


It is apparent from Fig. 4 that there is a close similarity between the weighted 
distributions of radiants recorded at different times of the year. Their most 
important features are therefore preserved in the combined distributions of 
Figs. 5 and 6. 

Fig. 5 shows the radiants to be concentrated in ecliptic longitude into three 
distinct groups. One of these lies at the longitude of the apex, and the other 
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two on either side of it at longitudes distant about 60°. The last two groups 
are particularly prominent at low ecliptic latitudes. 

The weighted latitude distribution (Fig. 6) is nearly symmetrical about 
the ecliptic in the observable range of /. An ecliptic concentration of radiants 
is apparent between latitudes 30°S and 30°N, and 50 per cent of the radiants 
in the weighted distribution lic within 15° of the ecliptic. North of latitude 
30° the weighted number of radiants per unit area is sensibly constant. 


Fic. §.--Distribution of radiants in ecliptic longitude. Meteors with radiants north of ecliptic only, 
(a) observed distribution 
(6) distribution corrected for observational selection. 
Curve 1, ecliptic latitude less than 10° 
Curve 2, ecliptic latitude less than 20° 
Curve 4, ecliptic latitude less than 40° 
Curve 4, ecliptic latitude less than go” 
Abscissa: longitude measured from the apex of the Earth's way. 
Ordinate: percentage of total distribution per 10° of longitude 


The radiant distribution, before inclusion of the weighting factors, is not 
asymmetrical about latitude o°, because the apparatus is mest sensitive to meteors 
whose radiants have declinations around 45°. ‘The symmetry of the weighted 
distributions may be considered evidence for the correctness of the weighting 
factor W,. 

(ii) Comparison with previous results.-A comparison of these results with 
those obtained by photographing bright meteors is unnecessary, because any 
differences arise from ditlerences between the orbits and these are discussed 
in Section 7. 

Comparison is of use, however, with the statistical data obtained by Hawkins 
(6), using a radio-echo apparatus only slightly less sensitive than the present 
one. Hawkins measured the radio-echo rates on two aerials directed at 
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different azimuths. He deduced that sporadic meteor radiants are concentrated 
towards the apex, and also towards two points on the ecliptic distant about 
60° from it. ‘The present radiant distributions exhibit these features, and are 
consistent with Hawkins’s results. 


Fic. 6.—-Distribution of radiants in ecliptic latitude. 
--- observed distribution. 
- distribution corrected for observational selection. 
Abscissa: ecliptic latitude. 
Ordinate: percentage of meteors with radiants north of the ecliptic per 5° of latitude, 


(b) The velocities, (i) The measured distribution.—The distribution of the 
measured velocities of the survey meteors is shown in Fig. 7. Weighted and 
unweighted distributions are given for all the meteors observed in the survey, 
and also for those with radiants north of the ecliptic. 

There is a prominent peak in the distribution around 30km/sec, a minimum 
between 45 and sokm/sec and a second smaller peak near 60km/sec. The 
number of meteors with velocities greater than 60km/sec is small and none is 
recorded with velocity exceeding 72 km/sec (the maximum velocity for a meteor 
belonging to the solar system). 

(ii) Comparison with previous results.—The velocity distributions obtained 
photographically by Whipple also show two peaks. ‘These occur near 25 and 
65 km/sec, and are more pronounced than those in the present results. ‘The 
gap between 45 and 50km/sec is also more prominent, and meteors having 
velocities less than 20km/sec, or greater than 60km/sec, are relatively more 


| 
10 
Pde 3 
4 ‘ 
; 
A : 
% : 
‘ 
i 
\ 
#0 4o o 60 


452 Jj. G. Davies and Fj. C. Gill Vol. 121 


numerous. It is shown in Section 7 that the differences between the velocity 
distributions are due te differences between the orbits of bright and faint 
meteors. 
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Fic. 7.—Distribution of observed velocities. 
(a) observed distribution. 
(6) distribution corrected for observational selection. 
Curve (1) all meteors. 
Curve (2) meteors with radiants north of the ecliptic. 
Abscissae: velocity in km|sec. 
Ordinates: percentage of meteors with radiants north of the ecliptic per 2 km|sec velocity 
interval. 

The velocities measured in the survey may also be compared with the 
radio-echo results of McKinley (5), and of Almond, Davies and Lovell (4). 
McKinley’s apparatus viewed approximately the same region of sky as the 
present one, and had much the same sensitivity. When allowance is made for 
the effects of shower meteors, and of instrumental selection favouring observation 
of the faster meteors, it is found that there is no significant difference between 
his results and those here obtained. Almond, Davies and Lovell observed the 


meteors coming from near the apex, with the object of finding out whether or 
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not any had velocity exceeding 72km/sec. In the autumn of 1951 they used 
the same transmitter and transmitting aerial as the present apparatus, and 
therefore viewed almost the same region of the sky. Their velocity distributions 
are found not to differ significantly from those recorded with the three station 
apparatus at the same times of day and year. 

The results of the survey thus do not conflict with previous observations of 
meteors of similar magnitude. 


6. The measured orbits 

The survey provided orbital measurements of 2474 meteors. Of these, 
the 1997 that had radiants in the hemisphere north of the ecliptic are included 
in the analysis below. 

The parameters computed for each orbit included: 


q, the perihelion distance, in astronomical units. 
q’, the aphelion distance. 

a, the semi-major axis; and 1/a. 

e, the orbital eccentricity. 

i, the inclination of the orbit to the ecliptic plane. 
w, the argument of perihelion. 

Q the longitude of the ascending node. 


The value of 1/a is a convenient measure of the energy of the orbit; g, q’ 
and e give its size and shape. The orbital plane is defined by the angles i and 
i, the latter angle being determined by the Earth’s position at the time of 
observation. The distributions of 1/a, e and i for the measured orbits are given 
in Figs. 8,gand1o. In Fig. 11 is shown the distribution of #, classified according 
to a and e. The most important features of the distributions are summarized 
below. 


(a) The distribution of 1/a.—It is seen from Fig. 8 that for almost all the 
orbits 1/a is positive; the orbits are therefore elliptical. ‘The inclusion of the 
weighting factors has little effect on the distribution. For about 2 per cent 
of the meteors, the calculated orbits are hyperbolic (1/a negative); but in 
every case the negative value of 1/a can be ascribed to the effects on an elliptical 
orbit of the anticipated errors of measurement. 

Most of the orbits have 1/a much greater than zero. For 60 per cent it is 
greater than o-5 (period less than 2-8 years), and for 22 per cent it exceeds 
10 (i.e. period less than a year). 


(b) Distribution of eccentricity.—The eccentricities of the orbits are shown 
in Fig. 9. Elliptical orbits of high eccentricity are the most abundant, but a 
substantial number of the meteors moved in orbits which were nearly circular. 
The way in which the eccentricity is related to the orbital energy and inclination 
may be seen from Fig. 11. 


(c) Distribution of inclination This is given in Fig. 10. Most of the 
meteors intercepted by the Earth move in direct orbits (i less than go°); there 
is a sharp peak at low inclinations, and a significant lack of orbits with inclination 
around go”. 

When the astronomical selection factor is included, to give the true distribu- 
tion of meteors in those orbits which can be intercepted, it is found that the 
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concentration near i =o ° disappears, the predominance of direct orbits becomes 
greater, and the lack of orbits with i near go° remains. 
If the orbits were randomly distributed, the number at inclination 7 would 


be proportional to sin 7. Orbits with i near go° would therefore be the most 
common. 


Fic. 8.—Distribution of reciprocal semi-major axis. Meteors with radiants north of ecliptic. 
(a) observed distribution 
(6) distribution corrected for observational selection 
(c) distribution corrected for observational and astronomical selection. 
Abscissae: reciprocal semi-major axis, t/a. 
Ordinates: percentage of distribution per 0-1 interval of 1/a. 


(d) The distribution of 1/a, i and e.—In Fig. 11, the orbits have been divided 
into four ranges of 1/a, and for each of these the inclination distribution is given 
for several ranges of e. 

The longer period orbits (with a@>2a.u.) intercepted by the Earth are 
concentrated towards i=o°, and are scarce near i=go° and near i=180°. 
The ecliptic concentration is most pronounced for the most eccentric orbits. 
When the astronomical selection is allowed for, the proportion of meteors in 


orbits of inclination less than 20° is decreased; it is, however, still greater than 
if the distribution were random. 
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Fic. 9.—Distribution of eccentricity. Meteors with radiants north of the ecliptic. 
Circles: observed distribution. 
Squares: distribution corrected for observational selection. 
Crosses: distribution corrected for observational and astronomical selection. 
Abscissa: eccentricity. 
Ordinate: percentage of distribution per 0:1 interval of eccentricity. 


Fic. 10.—Distribution of inclination. Meteors with radiants north of the ecliptic. 
Circles: observational distribution. 
Squares: distribution corrected for observational selection. 
Crosses: distribution corrected for observational and astronomical selection. 
Abscissa: inclination, degrees. 
Ordinate: percentage of distribution per 10° of inclination. 
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The more eccentric of the orbits with semi-major axes less than 2a.u. are 
distributed in a similar manner. The orbits with eccentricity less than os, 
which are observable only when a is sufficiently small, are rarely found at low 
inclinations. Instead, they are most abundant around inclinations 60° and 
140°, the grouping being most apparent amongst the orbits of shortest period. 
Again, few orbits are found with 7 near go’. 
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Fic. 11a.—Distribution of orbits as a function of semi-major axis, eccentricity and inclination. 
Observed distribution; meteors with radiants north of the ecliptic. 

Curve (1) e€<1°0 
(2) e<o'7 
(3) e<ors 
(4) e<or2. 

Abscissae: inclination, degrees. 

Ordinates: percentage of distribution per 10° of inclination. 


7. Discussion 

(a) The absence of interstellar meteors.—Of the 2474 meteors observed in this 
survey none was established to move in a hyperbolic orbit. It is thus unlikely 
that more than 0-05 per cent of observed meteors come from outside the solar 
system. 

(b) Meteor showers detected in the survey.—The 6-Aquarid shower, from 
which 35 meteors were recorded in 1954 July 27-28, was the only major shower 
observed in the survey. The radiant of this shower is south of the ecliptic, and 
therefore does not affect the distributions given above. Although no 
exhaustive search has been made for statistically significant groupings of 
meteors, it is clear that there are few other showers distinguishable from the 
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radiant and velocity data. Details of the showers that were detected are given 
in Table IV. 


TABLE IV 
Showers detected during the survey 

Number and Name I a 3 

Daytime 

Arietids Aquarids 
Radiant R.A. (deg.) 340 50 339 119 
Radiant Dec. (deg.) —5§ +26 —19 +15 
Observed vel. (km/sec) 29 4! 43 4I 
Perihelion distance (a.u.) 
Eccentricity 
Inclination (deg.) 
Arg. of perihelion (deg.) 
Long. of node (deg.) 
(deg.) 
No. of meteors observed 


Of the six showers listed, three (numbers 1, 4 and 6) were not previously 
known. Their orbits have low inclinations to the ecliptic and are not at all 
unusual. The similarity between orbits 4 and 6 is remarkable. It is also of 
interest to note that one meteor was recorded in 1955 March in an orbit closely 
resembling those of Whipple’s Virginid meteors (2). The orbits are compared 
in Table V. 


TABLE V 
Comparison of meteor 2537 with photographic Virginids 
Radio meteor 2537 Harvard meteors 
2905 1934 828 1158 

Radiant R.A. (deg.) 186:0 175°6 189°7 183°0 186°5 
Radiant Dec. (deg.) +3°0 +7°1 
Heliocentric Vel. (km/sec) 37°8 42°1 
Perihelion distance (a.u.) 0°38 
Eccentricity 0°85 
Arg. of perihelion (deg.) 289°8 
Long. of node (deg.) 38772 
a (deg.) , 287°0 270°7 277°3 


It has been noticed by several observers (e.g Kaiser (15)) that meteor showers 
are much less prominent amongst the faint than amongst the bright meteors. 
This is borne out by the fact that the four showers in Table IV (Nos. 1, 2, 5 
and 6) that had radiants north of the ecliptic contributed only 31 meteors, or 
about 1} per cent of the total number, towards the orbit analysis. Although 
some additional weak showers may be present, it is clear that very few of the - 
observed meteors belonged to showers, and that no allowance need be made 
in the analysis for their presence. 

(c) Comparison with the photographic results—The 51 orbits of sporadic 
meteors (all brighter than zero magnitude) that have been published by Whipple 
(2) fall into two distinct groups. In the larger group, estimated to contain 
about 80 per cent of the meteors intercepted by the Earth, the orbits are similar 
to those of the short-period comets. They are direct-moving, with inclinations 
usually less than 20°, and their aphelia lie near the orbit of Jupiter. The 
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orbits of the remaining 20 per cent of the meteors are almost parabolic and, as 
far as can be told from a small number, are randomly distributed in inclination. 
They are thus similar to the long-period cometary orbits. The data are too 
few to decide whether either group lacks meteors with orbits inclined at go°. 

In the present survey the orbits with semi-major axes between 2 and 5 a.u. 
are, like Whipple’s, mainly direct-moving and concentrated towards low inclina- 
tions. No distinct group of near-parabolic orbits is distinguishable amongst 
the radio results. This may be to some extent a consequence of the inferior 
accuracy of the radio method of measurement. 
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Fic. 11),—Distribution of orbits as a function of semi-major axis, eccentricity and inclination. 
Distribution corrected for observational selection; meteors with radiants north of the 
ecliptic. 
Curve (1) e< 1-0 
(2) e<o-7 
(3) e<ors 
(4) e<o-2. 
Abscissae: inclination, degrees, 
Ordinates: percentage of distribution per 10° of inclination. 


The radio results differ significantly from those of Whipple, however, in 
including many orbits of very short period (i.e. with semi-major axes less than 
2 a.u.). These are absent from the photographic results. 
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The possibility that the short-period radio orbits arose from systematic 
errors (of 10 or 20km/sec) in the measurement of velocity can be discounted, 
in view of the agreement with previous results of the measured Geminid and 
5-Aquarid orbits. Processes of selection, which could conceivably favour 
observation of short-period meteors, were also considered. Such selection 
would happen if, for example, these meteors were less likelythan others to break 
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Fic. 11¢.—Distribution of orbits as a function of semi-major axis, eccentricity and inclination, 
Distribution corrected for observational and astronomical selection; meteors with 
radiants north of ecliptic. 
Curve (1) e<1°0 
(2) e<o-7 
(3) e<os 
(4) 
Abscissae: inclination, degrees, 
Ordinates: percentage of distribution per 10° of inclination. 


into fragments on entering the atmosphere. Approximate orbit measurements 
(similar to the velocity measurements of section 3 (c)) were therefore made for 
80 per cent of the meteors recorded during the period 1954 May 20, 1547-2352 
sidereal time. The distribution of these orbits (approximately 400 in number) 
was not significantly different from the weighted distribution of those, amongst 
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them, that were included in the survey. Thus no evidence for a selection 
process was found. | 

(d) Radio-echo measurements of bright meteors.—The differences between 
the radio-echo and photographic surveys can be attributed to the difference 
between the magnitudes of the observed meteors. 

Further evidence that the orbital distribution of sporadic meteors depends 
on their magnitude is given in Fig. 12. This shows the distribution of 1/a 
amongst 55 meteors which were recorded during the survey, but were excluded 
from the main analysis because they gave ‘long-duration’ echoes. These 
echoes come from trails with more than 2 x 10! electrons per cm, and the meteors 
producing them were at least three magnitudes brighter than those in the main 
survey. Their orbits nearly all have 1/a less than 0-5, and are very like those 
of Whipple’s bright meteors. 
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Fic. 12.— Distribution of reciprocal semi-major axis for 55 meteors of magnitude brighter tian + 4. 
Abscissa: reciprocal semi-major axis, t/a. 
Ordinate: percentage of distribution per 0-1 interval of 1/a. 


(e) A possible explanation of the results of the survey.—A qualitative explana- 
tion of the survey results, and of the differences between them and those for 
bright meteors, is given below in terms of the action on meteors in long-period 
orbits of solar radiation and capture and perturbation by the planets. 

Consider the effect of solar radiation on meteors initially in long-period orbits 
that are randomly distributed in inclination and most of which lie wholly outside — 
Jupiter’s orbit. 

The effect of radiation on the orbit of a small body has been discussed by 
Poynting (16) and by Robertson (17). It is found that a body in an orbit 
round the Sun gradually loses angular momentum so that, if it is large enough 
not to be driven away by radiation pressure, it eventually falls into the Sun. 
Because the momentum is lost mainly when the particle is near perihelion, only 
the aphelion distance at first decreases. When, as a result of this, the orbit 
has become almost circular, it decreases rapidly in size and the particle spirals 
into the Sun. For a given orbit the rate of drift is proportional to the ratio of 
the surface area of the particle to its mass. 

A stage is therefore reached by a meteor originally in an orbit outside Jupiter’s 
when the orbit has become nearly circular and close approaches to the planet 
are possible. 

Opik (18) has discussed the probability of collision of a meteor, in a con- 
tracting orbit, with one of the planets. The probability is much larger for Jupiter 
than for any other planet, and is greatest for meteors in orbits that lie near the 
plane of Jupiter’s orbit (approximately the ecliptic plane). It increases with 
the size of the meteor. 
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An extrapolation of Opik’s results, which refer to meteor orbits of low 
eccentricity and inclination 10°, has been made to include other inclinations. 
If the surface area/mass ratios for the radio and photographic meteors may be 
assumed to be 50cm?/gm and 5 cm?/gm respectively, the destruction of almost 
all the photographic meteors with orbital inclination less than 40°, or greater 
than 165° will be likely by collision with Jupiter. For radio meteors, which 
drift more rapidly towards the Sun, the figures are 15° and 175 °. 

Apart from those actually captured, many more meteors will be 
perturbed into modified orbits. Faint meteors in orbits of inclination 
near go° have the smallest chance of collision, and may often reach orbits 
wholly within Jupiter’s after suffering only minor perturbations. The pertur- 
bations will tend to reduce the angle between the orbital plane and the ecliptic, 
so that the number of meteors in orbits with inclination go” will be diminished. 
When their orbits have contracted further, so that collisions with the earth are 
possible, the meteors may be observed. The orbits will be nearly circular, and 
have inclinations well away from the ecliptic. 

Bright meteors, and faint meteors in orbits lying closer to the plane of the 
ecliptic, are likely to suffer major perturbation before their orbits can contract 
sufficiently to escape from Jupiter’s vicinity. The perturbed orbits may be 
expected to be predominantly direct and of low inclination, and to have aphelia 
near Jupiter’s orbit, in the same way as do the orbits of short-period comets. 
If the orbits are sufficiently eccentric, so that their perihelia are within the 
Earth’s orbit, they may be observed without further change. 

In the case of the bright meteors the Poynting-Robertson drift is sufficiently 
slow for a ‘Jupiter family’ of orbits, similar to those of comets, to be formed. 
The planet then retards the further contraction of the orbits, so that at any time 
there are likely to be few bright meteors in orbits of shorter period. Faint 
meteors also may be retained in the group as long as their perihelia are far from 
the Sun so that their drift is slow. That these orbits are not observed 
probably indicates that the perihelion has to be outside the Earth’s orbit for 
this to happen. Faint meteors in orbits with short perihelion distances will not 
be retained by the group, and may be observed at any stage of their subsequent 
contraction until their aphelia pass within the orbit of the Earth. 

On the basis of this model, the following orbits might thus be expected: 


For bright meteors 

(i) Some randomly distributed orbits of long period, which have not yet 
suffered close approaches to Jupiter. 

(ii) A group of orbits like those of the ‘short-period’ comets, which have 
been formed from the long-period orbits by the action of Jupiter. 


For faint meteors 

(i) Some long-period orbits. 

(ii) Orbits of very short period, that are almost circular. These should be 
rare at inclinations 0°, go° and 180°. 

(iii) Direct orbits of low inclination and high eccentricity, many with very 
short periods, formed from long-period orbits after major perturbation by 
Jupiter. The ‘Jupiter group’ of orbits will not be prominent. 

The main differences between the photographic and radio results thus 
appear to be explicable, at least qualitatively, in this way. 
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8. Conclusion 

The use of the radio-echo method to measure the orbits of faint meteors 
has been demonstrated, and it has been found that the distribution of sporadic 
meteor orbits is dependent on the size of the meteors. A possible explanation 
of the dependence has been advanced. 

It is intended to use the radio-echo method of orbit measurement to 
investigate the phenomenon more fully. Measurements of the orbits of meteor 
showers will also be made. 
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SPECTRA OF A BRIGHT CORONAL PROMINENCE 
ASSOCIATED WITH A CLASS 3 FLARE 


Ian Elliott, M. A. Ellison and John H. Reid 
(Received 1960 July 16) 


Summary 


Photographs and spectra are described of a prominence which formed in 
the corona following a Class 3 flare at the limb. Line profiles have been 
obtained which are Gaussian in shape. On the assumption that the line 
broadening is of purely thermal origin a kinetic temperature ~ 10° °K is 
indicated for the Ha line. ‘The line-widths of H and K require temperatures 
~1o® K. It is concluded that these great line-widths must arise largely 
from non-thermal processes occurring within the prominence. Some possible 
mechanisms are discussed. 


1. Introduction.—A brilliant arch prominence associated with a Class 3 flare 
near the limb was observed with the combined spectrohelioscope and spectro- 
graph of the Royal Observatory, Edinburgh, on 1957 April 16. Two spectrum 
plates of the prominence were obtained in the H and K and Ha regions and the 
spectrophotometric results are here considered in relation to those for normal 
quiescent prominences. 

2. Flare and prominence.—The flare was first observed at 10" 40" U.T. on 
the east limb in latitude + 31° over a large complex spot group (long. 260°; 
lat. + 28°) which was just appearing. This region (HAO Specific Region 57-QA) 
showed complex coronal loop prominences for two days at the limb and was a 
source of impressive flares during the early part of its disk passage. The spot 
group declined rapidly and was last seen on April 27, five days after CMP. 
Meudon reported the flare as beginning at 10" 40™, maximum intensity at 11"05™ 
and ending at 13"00™. Partial observations were made at eight other stations (1). 

Up till the time of maximum intensity the flare emission had the appearance 
of half-a-dozen bright nodules some of which were seen against the disk, others 
against the sky, but none exceeding some 10 x 10°km in height. By 11°28" 
an arch had formed over the flare emission regions, the top of the arch being 
at 20 x 10°km above the limb and connected to the chromosphere by thin bright 
threads at its north and south ends. The photographs of Plate I illustrate the 
development of the prominence in the interval 11517™ to 12"01™. Through 
the kindness of Professor R. Coutrez and Dr A. Koeckelenbergh these prints 
were enlarged for us from the film taken with the Ha heliograph at Uccle. It may 
be remarked that only limb prominences of exceptional brightness show up on 
normal exposures with the heliograph, the exposure time being controlled to 
give the best record of disk features. During the next hour the arch changed 
little in brightness but rose slowly to 50 x 10°km (mean velocity across the line 
of sight 8km/sec). At this time (12" 30") it appeared as five or six closely bunched 
loops, extending over 7° of latitude, the southern ends of which terminated in 


? 
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the chromosphere close to the spot group (see our drawings in Plate II). by 
13" 25™ the arches had faded and were no brighter than the material of a normal 
quiescent prominence. At 14% 40™ a faint portion was visible at a height of 
70 x 10°km and at 16" 33™ very little of the arch structure remained. 

3. Terrestrial effects —Sudden ionospheric disturbances typical of a Class 3 
flare were reported from many sources. A sudden enhancement of atmospherics 
(S.E.A.), recorded at Edinburgh on a frequency of 24 Kc/s, began at 10" 40", 
reached maximum at 10" 50™ and ended at 11"47™. The receiving stations of 
Cable and Wireless Limited reported 54 short wave channels affected, with 
complete fade-outs (S.W.F.) in go per cent of these. A sudden cosmic noise 
anomaly (S.C.A.) was recorded at Edinburgh on frequency of 18 Mc/s with 
maximum absorption at 10"54". Geomagnetic crochets were reported from 
many magnetic observatories with maximum amplitude occurring at about 
10°50". Major bursts of solar radio emission were registered on 3000 Mc/s 
(10 40"—11" 30™), on 500 Mc/s (10" 46™—11" 02™), and on 200 Mc/s (10" 47""— 
11" 03") (2). 

We wish to emphasize that these terrestrial effects were associated with the 
flare emission which was observed at a low level rather than with the prominence, 
as they had either passed their maximum, or had returned to normal before the 
prominence developed. The emission of the prominence was, therefore, 
qualitatively different from that of the flare and it seems to have made no contri- 
bution to the ionizing radiations or to the radio noise emissions. 

4. Observational material.—'Two spectral plates were obtained at 11" 48™ (1) 
and 12" 11™ (2) with the spectrograph slit crossing the bright top of the prominence 
arch at 35 x 10°km above the limb. Plate No. 1, taken on Kodak Oa-O emulsion, 
includes the lines H8, He, H&, H, K and Het 4026.3A. Plate No. 2, on Kodak 
I1I-F emulsion, records Hx and He1 6678 A. 

Immediately after each exposure the plates were calibrated with the same 
exposure times as for the prominence (No. 1, 30 secs.; No. 2, 40 secs.) in the 
spectral light from the centre of the solar disk, using a platinum sputtered step 
weakener plus neutral screens of known transmissions for the respective wave- 
lengths. ‘This procedure enables us to express all line intensities as a fraction 
of the true continuum of the disk centre, which we call absolute intensities. 

5. Photometry.—The plates were brush developed in metolhydroquinone 
for 5 minutes at 18°C. The principal emission lines are shown in the enlargements 
of Plate II. 

The plates were measured with the Kipp and Zonen recording microphoto- 
meter of Dunsink Observatory. Two traverses were made across each plate, the 
first to record the emission line profiles at points A, A’, B, C, D (Plate 11) and the 
second at a point outside the prominence image to record the background of 
scattered light from the sky. By keeping the ‘‘lead-screw ’’ locked to the paper 
drum during these two operations both tracings were superimposed in wave- 
length register on the one sheet of bromide paper. 

The linear dispersions on the plates were: Plate 1 (1st order) 3:11 Ajmm 
and Plate 2 (2nd order) 1-47 A/mm. The length of the analysing slit of the 
microphotometer corresponded to a distance of 6000 km on the solar disk, and 
the analysing slit width (15) was equivalent to 0-047 A for Plate 1 and 0-022 A 
for Plate 2. 
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PLate I.—Development of the coronal prominence following a limb flare 
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The calibration spectrum on each plate, referring to the centre of the solar 
disk, was traversed on both sides of the emission lines at wave-lengths that were 
effectively free from Fraunhofer absorption lines. These wave-lengths and the 
effective continuum intensities are given in a previous publication (4, Table 4). 
In the region of the H and K lines the choice of a suitable continuum was rendered 
difficult by the presence of faint absorption lines and it was necessary here to 
estimate an ‘‘ effective continuum intensity ’’ by reference to the Utrecht Photo- 
metric Atlas of the Solar Spectrum. 

The conversion of photographic densities into light intensities was carried 
out by the methods which have previously been described in the Edinburgh solar 
spectrophotometric surveys (3, 4). Finally, the measured intensities of the 
emission lines have been corrected for grating ghosts (3). Corrections for 
instrumental profile have been applied to the observed profiles H« (A) and K (D) 
(see Fig. 2). 

6. Equivalent widths and central intensities compared with those of quiescent 
prominences.—The measured profiles of the emission lines and of the backgrounds 
of scattered light were first plotted, as described in a previous publication (4, 
p. 38 and Fig. 4). The corrected profiles shown in Fig. 1 were formed by sub- 
tracting the backgrounds in each case, and the equivalent widths (W’A) were 
determined by measuring the areas under the profiles with an Amsler planimeter. 
In Table I we give the equivalent width (W’A), the half-width and central 
intensity for each line. For comparison, Table I also contains the mean values 
of the same quantities for a number of quiescent prominences, as determined 
in one of the Edinburgh surveys (4, Tables 5 and 6), the same instrument and 
methods of reduction having been used in both cases. 


TABLE I 


Measured central intensity. 

Unit=10~* x continuum of 

centre of disk at same 
wavelength 


; Line width in A 
Plate Equivalent width 
No. Spectral line WA in central 


intensity 


(A) 122 
Hod (A’) 528 2°55 
He 1 6678 (B) 26 0°88 


H6(C) o'72 
He (C) 0°74 
H(C) 0°68 
K(D) 


Quiescent prominences (mean values) 


0°73 
0°43 
0°49 


The equivalent widths of H and K in the coronal prominence are ~ 35 and 
the central intensities ~20 times those of the quiescent prominences. These 
values seem to be much in excess of any previously recorded. The central 
intensity of the K line, for example, is 27 per cent of the intensity of the adjacent 


105 
2 
19 
45 
52 
222. 
274 
No. of 
cases 
13 Hx 44°6 57°8 
7 18 41 
3 He 16 3°5 
8 H are) 
8 K 
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Fic. 1.—Plotted profiles of emission lines. 


Ordinates: Intensities in terms of the continuum at the same wave-length for the centre of the disk. 
Abscissae: Angstroms. 
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continuum as measured at the centre of the disk. The H« region was photographed 
23 minutes later, at a time when the prominence had considerably faded: the 
equivalent width of Ha (A’) was then 12 times and the central intensity 3-5 © 
times the mean values for the quiescent prominences. We may note that the 
traverse of the Ha line at A gives a value for the central intensity (C.I.) which 
falls well on the plot of C.I. against W’A which was obtained previously for 
31 quiescent prominences (4, p. 43, Fig. 5). But for A’ the central intensity 
lies far beyond the extremity of this plot. ; 

7. Temperature estimates from the profiles.—The profile Ha(A) is Gaussian 
in shape and is closely fitted by a Maxwellian distribution of atomic velocities 
corresponding to T=50000°K (Fig. 2a). There is no appreciable flattening 
of this profile due to self-absorption, and so we assume an optically thin layer 
(i.e. optical depth C =o) in this part of the prominence. The slit was here crossing 
a down-flowing streamer and the profile has been Doppler shifted bodily to shorter 
wave-lengths (see Fig. 1) by 0-37A (sightline component of 17 km/sec. approach), 


10 K (D) 
T = 1,000,000° K 
C-5 


Fic. 2.—Prominence profiles compared with idealized temperature profiles. 


Ordinates: Intensities in terms of the central intensity as unity. 
Abscissae: Angstroms. 

(a) H(A). The observed profile of Fig. 1 has been corrected for instrumental broadening 
(crosses) and is compared with the profile for radiation from an optically thin layer (C=o) at 
T= 50 000 °K (full line). 

(6) Ha(A’). No correction has been applied for instrumental broadening owing to the great 
line-width. The observed profile (crosses) is compared with the temperature profile for 
T=150 000 °K and C=3 (full line). 

(c) K(D). The observed profile for the K line (Fig. 1) has been corrected for instrumental 
broadening (crosses) and is compared with the temperature profile for T=10°°K and C=5 
(full line). 


Notwithstanding the Doppler shift, this profile, like all the others, is highly 
symmetrical in shape. The profile Ha (A’), which corresponds in position with 
the extremely bright top of the arch, gives a fairly close fit to the temperature 
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profile of T= 150000°K, with optical depth C=3 (Fig. 2b). In other words, 
this profile is considerably flattened, possibly by self-absorption. However, the 
central density of the image here is rather too great for accurate photometry 
and we may have errors of the order of 10 per cent near the line centre. 

The H and K lines, recorded 23 minutes before Hx, are extremely well 
exposed and suitable for accurate photometry. We have corrected the observed 
K line profile for instrumental broadening, which has the effect of increasing its 
central intensity from 27 to 30 per cent. of the adjacent continuum for the centre 
of the disk. The resulting profile (Fig. 2c) is highly flattened and is reasonably 
well fitted except in the wings by a temperature profile for 7 = 10° °K and optical 
depth C=5. This state of affairs was quite unexpected and clearly underlines 
the necessity for some alternative explanation of the line broadening apart from 
purely thermal motions. 

8. Discussion.—The kinetic temperatures derived from hydrogen are un- 
doubtedly high, but they are of the same order as those obtained (100000°K) 
by Billings (5) at Climax for faint coronal condensation prominences occurring 
in highly active regions. Zirin (6), working with the 16-inch coronagraph and a 
high dispersion spectrograph at the U.A.R. Observatory, Sunspot, finds tempera- 
tures up to 50000°K in similar prominences. 

These ‘‘hot’’ prominences are of rare occurrence and the present example 
is perhaps exceptional because of its great brilliance. Kinetic temperatures are 
likely to be high in the coronal condensation regions, for the atoms here are in a 
state of transition from the extreme physical conditions in the corona (10° °K) 
to those within the stable quiescent prominences which consistently yield kinetic 
temperatures in the range 10000-20000°K (4). 

In these regions internal turbulence is probably high and possibly on a scale 
so small that it may be unresolved by the spectrograph. Under such conditions 
the line broadening may arise partly from thermal motions of the emitting atoms 
and partly from non-thermal motions of small turbulent elements, themselves 
having a Maxwellian distribution of sightline velocities. ‘The resulting line 
profile will still be Gaussian in shape, as we have found for Hx (A). 

However, something additional is needed if we are to account for the excep- 
tional line-widths observed in H and K. Here we are dealing with ions (Cal11) 
as opposed to neutral atoms (#/). It seems possible that the ions are spiralling 
along the lines of force of a magnetic field located within the prominence and 
they thus possess an added ‘‘ turbulence ’’ which the neutral atoms do not share. 
This kind of mechanism has recently been advanced by Athay (7) in order to 
explain the wide K profiles, having half-widths of o-8 A similar to ours, which 
he found in the spectra of chromospheric spicules. His kinetic temperatures 
derived from the lines D,, HB and Hy of the spicules were of the same order as 
those we have derived from Ha. 


g. Acknowledgments.—Our thanks are due to Professor R. Coutrez and 
Dr A. Koeckelenbergh for the use of their prominence photographs; and to 
the Astronomer Royal for Scotland for the observational facilities provided 
at Edinburgh. 
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RADIATIVE RECOMBINATION OF He 
A. Burgess and M. J. Seaton 
(Received 1960 August 4) 


Summary 


Recombination coefficients for 
He+ +e--He(n +hv 
are tabulated for T=o0, 1 x 10‘ and 2 x 10%. Hydrogenic data are used for L. >2, 


but allowance is made for the non-hydrogenic character of the He nS and nP 
states. 


Rate coefficients x(7) for radiative recombination may be calculated when 
the cross sections are known for the inverse process of photo-ionization (1). 
The coefticients for recombination of H*, 

H+ +e>H(nL) + hv (1) 
may be calculated exactly (2) and sums such as 
x(n) ='S and a(tot)= > a(n) (2) 
nel 
may be calculated, to an accuracy better than one per cent, using asymptotic 
expansions for the Gaunt factors (3). The «(#) are functions of temperature T 
and are such that 7"x(i) remains finite for 7->o but 7'*s(tot) diverges as log T. 
For recombination of He’, 
He+(1S) + e>He(n + hy, (3) 
the He states may be taken to be hydrogenic for L > 2, giving 
a(He,n aL) 
22, (4) 
a(He, = 3a(H, nL) 
but for 1 =o and 1 one must allow for the non-hydrogenic character of the He 
states. 
The photo-ionization cross section for the He1'S ground state is 
a,(He, 11S) = 7:3 x 107" exp ( — 0°76e) cm? (5) 
where the photon energy is 
hv =(1'807+€)/y, (6) 
I, being the hydrogen ionization potential and 1-807/,, the He ionization 
potential. For 


180714 Shy < (7) 

(5) gives agreement to within about one per cent with the best experimental 

and theoretical data (4,5). ‘Table | gives values of 7'*x(He, 1'S) obtained 

from (5) and, for comparison, values of T"°a(H,1S) (7 in deg. K and « in 
cm* sec '). 

For He P states and excited 5 states we have obtained photo-ionization 

cross sections from the general formula of Burgess and Seaton (4) and have 
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calculated recombination coefficients using the University of London Mercury 
Computer. Results are given in ‘Table I] together with comparison results 


for S and P states of H. ‘Table III gives results for the sums S a(n). 


TABLE I 


10% 7" "o(H, 1S) 


10!2 TJ" 2a(He, 1'S) 


16°46 15°6 
15°82 15°9 
15°26 16°2 


TABLE II 


Values of 10°T!?a 


H, nL He, n'L He, n°L 


A. A 


n L=o L=1 L=o L=1 L=o L=1 
2 2°413 6:634 *501 1°632 1:268 6°75 

3 ‘812 2°519 644 "295 2°25 

4 +381 1'206 ‘068 “112 1°020 
5 ‘213 ‘677 ‘031 ‘180 "055 "554 
6 ‘364 *422 ‘O19 ‘114 O31 339 
“ogo 284 ‘O12 ‘077 ‘020 "225 
8 ‘201 “008 ‘O55 013 158 
9 “149 ‘006 ‘O41 "009 “116 


‘005 


032 "007 


OC 


2x 104 


ANF 


| 

T=0 

10 037 114 089 

>> 4°41 12°80 “81 3°26 1°84 11°96 
nae 

T= 

1x 

2°342 5°355 "538 1°283 1°46 5°68 
-782 2°038 158 1°95 
. . . . 

‘ oe 359 966 076 246 153 886 
195 "532 ‘137 ‘079 475 
“118 324 ‘O22 084 047 285 
O77 ‘212 O14 054 ‘031 183 
"053 “146 “+010 ‘037 O21 125 
"105 ‘007 ‘027 ‘O15 “089 

‘078 ‘005 ‘020 ‘O12 065 
ie 4°13 10°13 ‘89 2°49 2°24 10°03 

T= 

2°262 4°588 "554 1068 1°57 4°93 
1°742 165 426 “414 1°72 
‘821 078 206 ‘172 
"449 "038 ‘114 "090 “419 
‘110 ‘272 "023 "069 "053 
‘071 ‘177 "044 "035 "159 
"049 ‘121 ‘O10 "024 "107 
‘035 ‘007 ‘021 ‘O17 ‘076 
"025 064 "005 ‘016 ‘013 

3°94 8°63 0°92 2°07 2°45 8°74 


Radiative recombination of He 


Tasce III 


Values of a(n) 
n=2 
H He, singlets He, triplets 
25°85 6°27 20°97 
20°19 4°89 16°91 


We wish to thank Mrs J. B. G. Wallace for her assistance with the computational 
work. 
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THE RELATIVISTIC MODEL OF THE STEADY-STATE 
UNIVERSE 


W. B. Bonnor 
(Received 1960 July 27) 


Summary 


It is shown that in the relativistic and Newtonian versions of the steady- 
state universe the motion of the cosmic matter is indeterminate. ‘The reason 
is that for matter with equation of state p +p =o both the inertial and passive 
gravitational mass densities vanish. ‘The conclusion is that these versions are 
of doubtful value in the prediction of cosmological observations. 


1. Introduction.—In the steady-state model of the universe space-time is 

described by a four-dimensional Riemann space with metric 

ds* = — + r? dé? +r? sin?6 dd?) + dt*, (1.1) 
where k is a constant, r, 6, ¢ are polar coordinates and tis the time. The specifica- 
tion of this metric does not by itself make the model complete. Field equations 
are needed to connect the expansion rate with the average cosmic density, and 
to work out the dynamics of individual galaxies. 

Attempts to provide a field theory have been made by Hoyle (1948, 1949, 
1960). Some advocates of the model do not accept Hoyle’s theory (Bondi 
1952), and in detailed investigations of the model it has not been customary to use 
Hoyle’s field equations (for example, Sciama 1955). In such work the tendency 
has been to use general relativity, or even Newtonian mechanics, and to suppose 
that these can be made to apply to the steady-state model. 

The most promising attempt to connect general relativity and the steady-state 
theory was that of McCrea (1951), and it is with McCrea’s work that this paper 
is chiefly concerned. McCrea’s theory has recently been the subject of an 
interesting investigation by Davidson (1959) and several of the results of this 
paper are similar to his. The conclusions drawn from the work described here 
are, however, quite different from those of Davidson. 

McCrea (1951) begins by noting that the field equations of general relativity, 
without cosmological term but with energy tensor 7%, viz., 

Ri — 38 R= (1.2) 
do not restrict the Riemann space to which they are applied, so it must be possible 
to use them to describe the space (1.1). To adopt this procedure amounts to 
using (1.2) to define an energy tensor of matter which would, according to 
relativity, produce the space-time (1.1) McCrea then shows that the matter 
concerned may be taken as having uniform positive density and negative pressure, 
and he concludes that the steady-state model may be described by the field equations 
of general relativity if the matter in it is assumed to have the necessary physical 
properties. 
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This procedure, at once simple and ingenious, seems unfortunately to fail 
on a crucial point. As will be shown, it so happens that for matter with the 
equation of state required by McCrea, viz., 


where p and p are the proper pressure and proper density in relativistic units, 
the motion is indeterminate. ‘Vhis weakens very considerably the predictions 
which can be made from the model, and makes it of doubtful value as an instrument 
of cosmological theory. 

This result is established in Section 2 and Section 3, and is followed in Section 4 
by a discussion of McCrea’s model in another coordinate system. In Section 5 
a Newtonian explanation of the indeterniinacy is attempted, and the paper 
ends with Section 6 in which the effect of the results on the steady-state theory 
is discussed. 

2. Motion of matter in McCrea’s model.—We assume that cosmic space-time: 
is described by (1.1) and use (1.2) to find the energy tensor of the matter which it 
contains. ‘The result (McCrea, 1951) is 


Ti = T3= T3= 3/82 ; (2.1) 


Let us form the eigenvalue equation at any event P corresponding to this 
energy tensor, viz. 
=0,(m=1, 2, 3,4). (2.2) 
where Y are the eigenvalues of 7} and w’ are the eigenvectors ; in this equation 
there is no summation over m. It is easily seen that the eigenvalues are all equal : 


{ 


(2) (3) 
X=X=X=X = 
Moreover, every vector is an eigenvector for each of the eigenvalues. Now 
the time-like unit eigenvector of an energy tensor at P gives the unit four-velocity 
vector of the matter present at P (Lichnerowicz 1955 ; Synge 1956). Since 
this eigenvector is arbitrary for 7* given by (2.1), it follows that the velocity of 
matter in McCrea’s model is completely undetermined. 

The usual procedure in the analysis of an energy tensor at a point is next to 

identify X, the eigenvalue corresponding to the time-like eigenvector, with 
(2) (3) 
the proper density, and —A, —X—NX with the principal pressures. In our 
case, this gives 
(dp (2) 


; 


The material of the model is then a perfect fluid with negative isotropic pressure 
p and density p such that 


p+p=o. (2.3) 
If we write the energy tensor for a perfect fluid in the usual form 
dx* dx' _ 
= Lia (P +p)—- dip, (2.4) 


we see that (2.3) allows the velocity vector dx'/ds to be arbitrary, as already found. 
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One, could of course assert arbitrarily that the matter in the model follows 
a geodesic, in which case one finds that its velocity is 


(2.5) 


and it is at rest in the coordinate system of (1.1). Reasons will be given in 
Sections 4 and 5 why this assumption would be unsatisfactory, but it may be 
noted here that one of the most powerful features of relativistic dynamics is 
that the field equations determine, in general, the motion of the matter present, 
and it would seem to be quite contrary to the spirit of the theory to assume 
here a motion which is not demanded by the tield equations. 

3. Significance of the indeterminacy of the motion.—TYo understand the 
significance of the statement that McCrea’s model does not determine the 
velocity of the matter present, we compare it with those of ordinary relativistic 
cosmology. The metrics for these are 

ds? = — + dd +r? sin*d dg*)(1 + far?) (3-1) 
and to be definite we suppose that one particular model is chosen so that g(t) 
and & are given. ‘Then, provided p+p#o, the field equations (1.2) together 
with the energy tensor (2.4) determine the velocity of the matter present to be 


and p and p are given in terms of g and its derivatives. ‘lhis shows that (3.1) 
refers to a coordinate system in which the particles in the model remain at rest. 
The model may now be used to make predictions about the red-shift, and other 
observable data. 

If this procedure is followed through for MecCrea’s model, then, as shown 
in Section 2, the velocity is not determined, and it is not known how particles 
are moving relative to the coordinate system being used (1.1). 

This relative motion is to a certain extent determinable by direct observation. 
For example, from the observed spherical symmetry one can assert that d@/ds 
and dd ds are zero. However, in the red-shift formula (and other similar 
formulae) there will now appear the unknown function dr/ds, depending on 
rand ¢t. ‘The effect of this is drastically to weaken the model as a means of 
prediction, since by appropriate choice of this function, it can accommodate 
widely ditferent observational possibilities. In fact, only a detailed investigation 
could show whether the model is capable of falsitication at all. 

4. The static form of the model.—The metric (1.1) is that of the de Sitter 
universe. It is known (Lemaitre 1925) that if new variables 7, / are introduced 
by 


(1 — 
(1.1) is transformed into 
ds? = — (1 — dv? — + sin?0 dd*) + (1 — (4.2) 


In this form, the space-time of the steady-state theory is static. 

If the field equations (1.2) are applied to the metric (4.2) we find that 7’ 
has the same values as in (2.1), and, as before, that the motion of the fluid is 
indeterminate. 


dr dO dt_ 
dr _ dd _db_ dt_, 
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We may, if we wish, assume that the matter is at rest in the coordinate system 

of (4.2), and take 
ds 

and this might seem to be reasonable because the model (4.2) is static. However, 
it is easy to see that if we adopt (2.5) for the model (1.1), we may not take 
(4.3) for the model (4.2). The reason is that (2.5) are the components of a 
four-vector, and when this is transformed according to (4.1) it does not lead to 
(4.3). Thus if we adopt McCrea’s assumption (2.5), the matter in the static 
model (4.2) must be moving with respect to the (7, #) coordinate system. In fact, 
of course, since (2.5) represents a geodesic, its transform must be a geodesic in 
(4.2). ‘This conclusion, that the matter in the static model (4.2) must be moving, 
is another argument against the gratuitous assumption that the matter follows a 
geodesic. 

Let us now turn to another curious circumstance, already noted by Davidson 
(1959). Consider the metric 


ds? = — (: — sin?6 d¢?* +(: =" dt?; (4.4) 
r 


this represents a universe similar to that described by (1.1) or (4.2), but with 
a mass m at the origin. Now, using (4.4) to evaluate (1.2) one finds again the 
same values for 7; as in (2.1), and the motion is arbitrary even in the presence 
of the mass. Jt seems that the matter with equation of state (2.3) is indifferent 
to a gravitational field. We shall return to this point in the next section. 

It should be emphasized that these strange conclusions apply to the model 
which contains this type of matter, and not to the de Sitter model. In the latter 
the field equations (1.2) are altered by the addition of the cosmological constant, 
and the energy tensor, the pressure and the density are all zero. A test particle, 
if introduced, follows a geodesic and its motion is in no way indeterminate. 

5. The physical interpretation.—To interpret physically the behaviour of 
matter in the model we start from the clue given towards the end of the last 
section, that matter satisfying (2.3) seems not to respond to gravitation. Let 
us adopt a Newtonian point of view, and consider the three ways in which mass 
density may appear in Newtonian mechanics. We define 

o, = inertial mass density ; 

0, = active gravitational mass density ; 

03 = passive gravitational mass density ; 
here active and passive gravitational mass refer to the power of a body to exert 
gravitation and to respond to a gravitational field in which it is placed. 

Of the three densities, only o, has an analogue in general relativity which is 
known with any assurance. In a static field it appears from the work of 
Whittaker (1935) that 

o,=T}-T}- 
and if matter is present in the form of a perfect fluid at rest in the coordinate 
system, we may write, with McCrea, 


Oe =p + 3p (5. I ) 
(relativistic units). 


| 
| 
| 
| 
is 
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To obtain a corresponding expression for o,, we note that for a perfect 
tluid without pressure 


dxi dx" 
= (5-2) 


and here p is clearly to be interpreted as the inertial mass density. Comparison 
of (2.4) and (5.2) suggests for a perfect fluid the tentative indentitication 


(5-3) 

Concerning o, very little is known. For a structureless point-particle 

it is a basic assumption of general relativity that the passive gravitational mass 

is equal to the inertial mass, but for extended material in the presence of stresses 

the equality by no means clearly follows. Nonetheless, | shall assume equality 
and write 


+p. (5-4) 
Evidence for (5.4) in the case of a perfect fluid at rest exists in the equation 
of hydrostatic equilibrium (Tolman 1934), which for the interior of a spherically 
symmetric body is 


dp d 
+ 2(p +P) log gas = 0, 


r being the radial coordinate. Taking } log g,, as the analogue of the Newtonian 
gravitational potential, we see that p + p is playing the part of the passive gravita- 
tional mass density. 

In this connection we may note the following interesting property of matter 
obeying the equation of state p+p=o0. Inserting this equation into the energy 
tensor (2.4) and forming the identities 

g=0 
we find, with no assumption whatever about the g,; except that they exist, 

¢ — 

so that the pressure (and density) must be constant throughout space-time. 
The impossibility of a pressure gradient is readily understood in the light of 
the preceding argument, since if the inertial mass density o, vanishes (as it 
will from (2.3) and (5.3)) a pressure gradient would produce an infinite 
acceleration. 

Let us now tentatively adopt (5.3) and (5.4). Then it follows from (2.3) 
that in McCrea’s model of the steady-state universe both the inertial and the passive 
gravitational mass densities are zero. For the Newtonian equation of motion of 
a small volume dv of matter in a gravitational field g, we have 

o, dof = 0, dug (5-5) 
f being the acceleration. Hence, if the densities o, and o; both vanish, f is 
indeterminate. We see that the indeterminacy of the motion has a quite natural 
interpretation on this extended Newtonian view. 

If this reasoning is correct, it follows that the Newtonian treatment of the 
steady-state model (McCrea 1951) is also doubtful. In forming the Newtonian 
equation of motion for the material substratum of the model, McCrea quite 


~ 
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naturally cancelled out what has here been called o, and 3. But if these are 
zero this is not permissible. In fact it seems that the motion of matter is also 
indeterminate in the Newtonian model. 

In general relativity the postulate that a test particle moves on a geodesic is 
equivalent to the Newtonian equation of motion in a gravitational field. To 
remove the indeterminacy in the relativistic model by postulating that the cosmic 
matter moves on a geodesic in (1.1) or (4.2) would amount, in the Newtonian 
analogue, to proscribing a particular fin (5.5). Since (5.5) evidently determines 
nothing if o, and o, are zero, 1 conclude once more that there is no justification 
for the postulate that matter moves on a geodesic. 

6. Conclusion. have argued in this paper that neither the Newtonian nor 
the relativistic versions of the steady-state theory, at least in their current forms, 
determine the motion of matter in the universe. Moreover, this lack of deter- 
mination weakens these versions to such an extent that they seem to be of little 
value as instruments of prediction in cosmology. 

It is just possible that the situation might be retrieved along the foliowing 
lines. ‘The indeterminacy in the motion arises only if 


pt+p=0 (6.1) 
is exactly satished. Could one therefore overcome the difficulty by assuming 
that (6.1) is only approximately realized in the actual universe? The answer 
is that the exact fulfilment of the equation is necessary if the universe is to be 
in a steady state, for the slightest departure from it will mean a secular change 
in the nature of the universe, which the authors of the theory are at pains to 
avoid. However, it might be possible to devise a model in which, although 
(6.1) is not satistied exactly in any large region of space-time, it 7s satisfied 
exactly over the model as a whole. Such a model would have different rates of 
expansion in different regions, and would be essentially inhomogeneous. It 
would probably be very complicated mathematically, and it might not accord 
with observation ; but the possibility of a model of this sort cannot be ruled 
out. 

The criticisms of the relativistic and Newtonian models given in the preceding 
sections do not, of course, apply to the steady-state theory itself. Since the 
latter does not accept general relativity or Newtonian mechanics on a 
cosmological scale, it is unaffected by defects in models constructed from these 
theories. What now seems clear, however, if my arguments are correct, is that 
those who wish to work with the steady-state theory must use a dynamics 
specifically designed for it, since they are not free to use the existing models of 
general relativity or Newtonian theory with any degree of confidence. 

To illustrate this point I may mention the theory of galaxy formation in 
the steady-state model. Interesting work with powerful conclusions has been 
done on this by Sciama (1955). However, Sciama uses Newtonian theory 
in a way which seems unjustified in the light of this paper. In particular the 
use of Newtonian mechanics to explain the condensation of newly-created 
cosmic matter in the presence of existing galaxies is not permissible if the 
inertial and passive gravitational masses of the created matter are zero. 
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Summary 


The intensities of the centre of the solar disk at 8-63 and 12°02 have been 
measured by a technique previously described for 11-104. The emission 
temperatures obtained were 5160 °K at 8°63 and 5050 °K at 12-02 against 
the earlier value of 5036 “K at 11-104. ‘The results are a few per cent greater 
than values calculated from model solar atmospheres. 

Limb darkening measurements have been made out to 0-91 of the solar 
radius (solar image diameter 17 mm) at 8-63, 11°10, 12°02 and 12°95, and are 
in good agreement with those of Pierce et al. at 8-3 and 10°24. 

The combined measurements of solar intensity and limb darkening have 
been used to study the wavelength variation of the photospheric opacity. 
The results show less variation with wavelength than predicted theoretically. 
They also indicate that the free-free component of H~ opacity calculated by 
Chandrasekhar and Breen is 40 per cent too high, which is approximately in 
agreement with the recent calculations of T. Ohmura and H. Ohmura. 


1. Introduction.—The position regarding solar photospheric opacity before 
the present work was undertaken was that the absorption coefficient of the 
solar atmosphere had been computed by Chandrasekhar and Breen (1) for bound- 
free and free—free electron transitions in the field of H~, the latter being dominant 
above 1-64. Neven and De Jager (2) and Peyturaux (3), using the above 
absorption coefficients together with a given model of the photosphere, have 
computed the central intensity and limb darkening for wavelengths up to 
10°24. The calculated central intensities and limb darkening up to 2-34 seem 
to agree fairly well with the observed values (3, 4, 5). However, beyond 2-34 
no central intensity measurements existed except that of Saiedy and Goody 
(6) at 11:10u. Limb darkening measurements at 8-34 and 10-2 were carried 
out by Pierce, McMath, Goldberg and Mohler (7) but their values are much 
less than those predicted, with Peyturaux’s model giving the better fit. 
Assuming the model adopted by Peyturaux (3) or Neven and De Jager (2) to 
be correct then the limb darkening measurements at 8-3 and 12-024 demand 
an increase in the absorption coefficients in this region, whereas more recent 
calculation of the free component of H~ absorption indicate the opposite (8). 

The present paper describes briefly new measurements of central intensity 
at 8-63 and 12°02u, and in more detail measurements of limb darkening at 
8-63, 11°10, 12°02 and 12-95, carried out at Ascot, Berkshire. ‘The combined 
measurements of central intensity and limb darkening provide a means of studying 
the spectral dependence of the photospheric opacity in the far infra-red. 
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Corrected air masses 
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Fic. 1a.—Results of observations on 8-63 
Filled circles: 1958 October 17 Open Circles: 1959 February 1 
Open squares: 1959 May 27 Filled squares: 1959 June 17 
Crosses: 1959 July 7. 
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Fic. 1b.— Results of observations on 12°02 
Filled circles: 1958 October 16 Open circles: 1959 April 28 
Open squares: 1959 May a Filled squares: 1959 May 5 
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2. Central intensity at 8-63 and 12-02u.—The measurements were carried 
out in the same way as those of an earlier experiment at 11-1ou reported in a 
recent paper (6). The window at 8-63. unfortunately contains N,O lines 
belonging to the 8-6. band and the absorption does not necessarily follow 
Lambert’s Law. However, the deviation from Lambert's Law at various air 
masses were obtained from unpublished data of the transmission of N,O at 
8-63 by Dr C. D. Walshaw in fitting an Elsasser model to laboratory measure- 
ments of Goody and Wormell (g). For these calculations the necessary 
parameters at 8°63 were «/5=0°188, where S=line intensity, 
a= Lorentz half-width and 5=line spacing, together with a constant volume 
mixing ratio of N,O to dry air of 3-5 x10~7 (10) and a mean temperature of 
250 K. 

The results of observations on the two windows, 8-63 and 12-02, are shown 
in Figs. 1a and 1b, and the computed central intensities, after applying a 
correction for limb darkening, in Table I with the estimated standard error. 
The equivalent central emission temperatures (7) were then calculated 
using Allen’s (11) data and are also included in Table I together with the results 
reported earlier. 


TABLE 


Central intensity at 8-63, 11°10 and 12.024 


Ip Standard error 


erg ster~! sec™! per cent 
per AA=1 cm 
6-530 x 0°87 
2-408 x 0°67 
1'774 


Tasie II 
Central emission temperatures from visible to the micro-wave region 
Authors 


I+ 


Peyturaux (3) 
Peyturaux (3) 
Peyturaux (3) 
Peyturaux (3) 
Peyturaux (3) 
Peyturaux (3) 
Present work 

Saiedy and Goody (6) 
Present work 

Coates (12) 

Coates (13) 
Whitehurst et al. (14) 
Coates (13) 
Zelinskaja et al. (15) 
Troickij et al. (16) 
Covington et al. (17) 


BS 


35 


A T> 
°K 
5036+ 30 
12°02 5050+ 80 
°K °K 
613 6190 
“883 6094 
1:008 6120 
1°546 6329 
2'029 6044 
2°310 5880 
8-63 5160 
5036 
12°02 5050 
4000 . 5600 
4300 6962 
6000 4500 
8600 8431 || 
3°2x 10! 13000 
10 x10! 50000 
10°3 x 10 21800 
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Table 11, summarizes the most up to date measurements of central emission 
temperatures from the visible to micro-wave regions, some of which are shown 
in Fig. 2. Fig. 2 shows a decrease in central emission temperatures and so of 
photospheric absorption coefficient (1) beyond 1-64. Around 5 mm there is 


0 t 


Fic. 2.—Variation of central emission temperature with wavelength. 


Open circles: — Peyturaux (3) Open squares: present work 
Open triangles: Whitehurst et al. (14) Filled circles: Coates (12, 13). 
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Fic. 3.—Observed and calculated central emission temperatures as function of wavelength. 
Filled circles: observation Open circles: Peyturaux (3) 
Crosses: Neven and De Jager (2). 
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a large uncertainty in the measurements of 7',._ However, our measurements 
fall very reasonably between those of Peyturaux in the near infra-red and that 
of Whitehurst, Coplan and Mitchel (14) at 6mm. Beyond 6mm there is a 
definite increase of 7. which indicates an increase in the chromospheric opacity. 

Our measurements of central intensity are compared with values calculated 
by Neven and De Jager and by Peyturaux at 8*3 and 1o-2y in Fig. 3. These 
authors used different photospheric models but the same absorption coefficient 
of H~ as calculated by Chandrasekhar and Breen. The inferred central 
emission temperatures 7. from the calculations of Neven and De Jager are on 
average 6 per cent lower than ours and decrease more rapidly with wavelength, 
while those of Peyturaux are systematically lower by 2 per cent. 


3. Limb darkening measurements between 8-6 and 13-0u 

3.1. Technique.-The observations were made with the spectrographic 
equipment (over-all time constant 0-7 secs) described by Saiedy and Goody 
(6). A siderostat, 30cm in diameter, and a paraboloid of equal diameter with a 
focal length of 180cm formed an image of the Sun, 17mm in diameter, on 
the entrance slit of the spectrometer. ‘The entrance slit (4:12 x0°25 mm) is 
vertical and consists of thick brass jaws in good thermal contact with the base of 
the spectrometer. This eliminates excessive heating of the slit, which could 
cause a variation in the slit size with time as the Sun drifts across it. The 
paraboloid was set so that its focus coincided with the entrance slit. This was 
achieved by moving the paraboloid backwards and forwards along the axis of 
the spectrometer until the image of a distant object (the Moon) coincided with 
the plane of the slit to within + 1mm at180cm. A distance-reading microscope 
was used in the final stages of this adjustment. 

Observations in the windows at 8-63, 11-10, 12°02 and 12°95 were carried 
out near noon on clear days without excessive convection, and the solar disk 
was scanned by stopping the siderostat drive. The vertical edge of the disk 
was set symmetrically, by eye, on the vertical edge of the slit using the declination 
and right ascension drives of the siderostat for fine adjustments. The image 
was then shifted horizontally off the slit and held there by the synchronous motor 
of the siderostat. After recording the zero level, the synchronous motor was 
switched off and the image drifted across the slit by the diurnal motion of the 
Earth. Deflections corresponding to the centre of the Sun were recorded 
before and after each drift record to allow for change in the sensitivity of the 
detector system. One such record taken at 8-63 is shown in Fig. 4. 

The scanning time was approximately two minutes. It is desirable to make 
the entrance slit as narrow as possible compared with the diameter of the image 
of the solar disk consistent with good signal to noise ratio. Since we do not 
aim at high spectral resolution, the exit slit may be quite wide, but in practice 
this was limited by the construction of the spectrometer and corresponds to 
approximately o-o5y at ron. The entrance slit width was such that a small 
feature of the solar image took two seconds to cross from one side to the other. 
The finite size of the entrance slit was taken into account by numerical integra- 
tions over the slit to obtain the effective value of cos@. Other factors that could 
distort the limb darkening curves are: 

(a) Non-linearity of the detector system, and stray radiation from shorter 


wavelengths. 
35* 


ge 
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(6) Lack of sharpness of the solar image. 

(c) The time constant of the detector and recorder. 

(d) Setting the solar image unsymmetrically on the slit. 

The effect of non-linearity and stray radiation from shorter wavelengths is 
completely negligible in our system. A diffuse solar image could arise from an 
out-of-focus slit, and by using the paraboloid off-axis, but focus was ensured 
to +1mm at 180cm and the effect of the off-axis paraboloid is to increase the 
height of the slit by approximately 4 mm, which could only be important on the 
very edge of the solar disk. ‘The time constant of the system was measured by 
applying an intensity step function to the detector and was found to be 0-7 sec. 
The asymmetry in setting the edge of the solar disk on the slit was up to } mm, 
for which reason the analysis of the records was restricted to cos@>0-4. In 
the region 1 > cos @> 0-4, the effect of the time constant is to shift the drift curve 
by 0-7 sec along the time axis (0-58 per cent of disk diameter) without changing 
its shape. ‘The records were measured at six points on each side of the centre 
of the drift curve and an average of each pair of the measurements was taken. 
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Fic. 4.—S is profile of limb darkening at 8°63. S,, Sy are deflections at the centre of the Sun; 
the decrease in signal at A,, Ay and Ay, is due to sun spots. The effect of the finite time constant 
could be clearly seen on S, and Sy. 

3-2. Results.—Five runs were made at 11-10 and 12-02, and four at 8-63 
and 12°95y, all in the spring of 1959. Table III gives the results, 


which are averages over all the runs at each wavelength, together with the 
standard errors. ‘There is a systematic increase of approximately 0-2 per cent 


Tasre III 
Results of observations on limb darkening. The listed values are 6(8) = 1(8)/I(O), 
= standard error representative of all observations 
12°95 
Cos (0) (0) 
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in limb darkening from 11-1 to 12-ou but no significance can be attached to this 
as it is of the order of the error in measuring the records. 

Fig. 5 shows our results in relation to those of Pierce et al. at 8-3 and 10-2, 
and reasonable agreement is found. 


7 8 9 10 i te 13 


Fic. 5.—Observed limb darkening in the far infra-red. 
Filled circles: present measurements Open circles: Pierce et al. (7). 


4. Opacity of the photosphere.—The angular and spectral distributions of 
radiation emerging from the photosphere give complete information about the 
radiation of each volume element within the photosphere (18). The emergent 
radiation at angle @ to the outward normal, assuming local thermal equlibrium, 
is given by 


1,(6)= | sec Odr, (1) 


where 
B,(T)= black body function for temperature 7 


7,= optical depth and is a function of T. 
When equation (1) is divided by the central intensity /,(©) then 


6,(0)= | bdr, (2) 
where 


BA(7,) 
S,(7,) = T(O) (3) 


For the study of the photospheric opacity S,(7,) is conventionally expressed 
analytically, so that equation (2) can be integrated. The simplest method which 
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fits our observations satisfactorily, originally proposed by Chalonge (19), is 
S,(7,)=4, (4) 


where a,, 6, and c, are constant at each wavelength. Then from equations 
(2) and (4) it can be shown that 
¢,(9) =a, +b, cos 6 + 2c, cos* 6. (5) 
A better representation has been given by Kourganoff (20) which has the 
correct asymptotic behaviour for both 7-0 and 7->o and is given by 
S,(7,)=@,' + (6) 
where 
E,(2)= | “<> du and 6,’ and c,' 
are also constant for each wavelength but differ from those given by equation 
(4). Integrating equations (2) with (6) gives 
$,(0)=a,' +6,’ cos +c,’ [1 —cos loge (1 + sec@)]. (7) 
Since our observations do not reach close to the limb there is no significant 
difference between (5) and (7) and the former has been used to provide, from 
Table III, by least squares, the values of the coefficients a,, 6, and c, given 
in Table IV. These lead, from (4), to the values of S,(7,) given in Table V. 


IV 


Coefficients a, b, and c from measurements of ¢,(0) 


a c 

“8865 — 
“9166 — "0405 
“9149 — 0382 
—°0344 


TABLE V 
Values of S, (7) from the observations of limb darkening and equation (4) 


A T 
8-63 “8865 92605 -9627 “9950 1°0479 1°OgS2 
“9166 *9478 ‘9758 1°0006 10404 10759 
12°02 "9149 9455 9730 *9975 1°0373 10740 
12°95 ‘9197 "9483 ‘9740 ‘9973 1°0352 10715 


4-1. Solar surface temperature T,.—TYhe solar surface temperature can be 

deduced from S,(7,=0) and J,(©) according to 
(8) 

and 

2hC (9) 
9 
but the values obtained depend very much on the form of S,(r,). Using 
S,(7,=0) from Table V and J,(©) from Table I, the surface temperatures 
T, given in Table VI under equations 4 and 5 are obtained. These values 
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should be considered with caution as they are an extrapolation of the experimental 
results. The variation with A is well within the error deriving from error in 
I,(©) and in a,, b, and c,. The surface temperatures inferred agree with those 
given by Peyturaux for the visible and near infra-red ; Peyturaux’s values ranged 
between 4356 and 4650°K. Our values also agree reasonably well with the 


VI 
Solar surface temperature T, 


A T, 
im °K eqns. (4), (5) °K eqns. (6), (7) 
8°63 4660 4332 
4667 4263 
12°02 4665 4379 


central emission temperature (4500 K) measured by Whitehurst et al. (14) 
for the mm range where the photospheric opacity results in emission only from 
the very outer layer of the photosphere. Kourganoff’s expression (eqn. 6) 
for S,(7,) has also been used to compute the surface temperature, the values of 
which are given in Table VI. ‘These show larger variability (but still within 
experimental error) than those from equations (4) and are less in the mean by 
340 °K. The decrease in surface temperature for this model is to be expected 
from the difference in form of (4) and (6) as r->o. 
4.2. Relation between average and central solar intensity.—TVhe average to 

central intensity is given by 


_ (x0) 
LO) 1,(O)fanrdr 
where /(r) =the solar intensity at a distance r from the centre of the solar disk. 
r=Rsin0@, R being the solar radius. 
The integration of (10) with (5) gives F,/1,(.0)=1—46,—c, and values of 
this ratio, using the coefficients of ‘Table IV, are given in Table VII. 


VII 
The ratio of mean to central solar intensity 
A F, 
1,(O) 

8-63 
12°02 0985 
12°95 0'986 


4.3. Determination of + as a function of 4.—The determination of the optical 
depth at different depths, detined by the temperature, enables us to compare the 
absorption coefficient for the photosphere at different wavelengths using 


B,(7,)=S,(7,) 1,(0) (11) 
with J,(©) from Table I and S,(7,) from Table V. The results of the computa- 
tion are shown in Fig. 6, the appropriate weight being given to each point in the 
curves drawn. It is evident that the absorption coefficient increases with wave- 
lengths as expected (1). However, it also appears that the increase is less than 


492 F. Saiedy Vol. 121 


predicted. For example, the experimental values of +7(11-10u)/7(8-63) and 
7(12°02,)/7(8°63) are approximately 1-36 and 1-38 for the whole upper layer of 
the photosphere (7 < 5300 °K) while the predicted ratios on the assumption that 
kad (1) in this region are 1-66 and 1-97 respectively. The errors in the above 
two observed ratios from the error in /,(©), around S,(r)= 1-0, were found to be 
13 per cent and 26 per cent, respectively showing that the spectral variation of the 
optical depth is just significantly less than the theoretical variation. ‘The large 
error in the second ratio reflects the uncertainty in /,(©) at 12-2, an error of 1 per 
cent in giving an error of approximately to per cent in r. 


4 


9 


Fic. 6.— Variation of the optical depth (+r) with wavelength (A). 


To compare our results with those obtained at shorter wavelength we have 
chosen the values of + measured by Peyturaux at 1-009 where Peyturaux’s and 
Mulders’ (21) values of central intensity agree best. The level, defined by 
T= 5300 K, was also chosen for this comparison for the following reasons: (i) 
Peyturaux's measurements of limb darkening extend only to cos @~ 0°16 which 
corresponds to r~ 0°16, i.e. 7'~ 5300°K and, for any lower value of T, + would 
strongly depend on the form of S,(7,) chosen; (ii) on the other hand, in the far 
infra-red 5300 K corresponds to + ~ 1-0 which is already a slight extrapolation of 
the experimental result. + at 5300 K was measured from Fig. 6 at 8-63, 11-10 and 
12°02 and expressed as a ratio (7, of the corresponding value of 1-oogu 
(Table V1 of Peyturaux). This ratio is compared with theoretical values of 
absorption coefficient in Fig. 7 in which Curve 1 derives from the bound-free 
and free-free absorption coefficient of Chandrasekhar and Breen, while Curve 2 
derives from the bound—free absorption coefficient of Chandrasekhar and Breen 
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and the free-free absorption coefficient of Ohmura and Ohmura. Our results 
agree better with the latter curve which indicates a smaller free-free absorption 
coefficient of H- than given earlier by Chandrasekhar and Breen. A better method 
of comparison between observations and theory is suggested by measuring the 
optical depth between two levels, say 5300 and 5000 K or 5300 and 5600 °K, in 
order to reduce the effect of any small error in /,(©) on the measured optical 
depth. This method is, however, unsatisfactory for the two reasons mentioned 
above. In this connection it would be a great help if measurements of central 
intensities around 4u were available. 


lt 


Fic. 7.—Experimental variation of the optical depth at 5300°K with wavelength, compared 
rotth theory. The value of + at 10094 was obtained from the measurements of Peyturaux. 
Curve 1: bound—free ~free—free (Chandrasekhar and Breen). 
Curve 2: bound—free (Chandrasekhar and Breen) free-free (Ohmura and Ohmura). 


5. Discussion.—Neven and De Jager and Peyturaux argued for absorption co- 
efficients larger than the theoretical values of Chandrasekhar and Breen in this 
region in conformity with the limb darkening measurements of Pierce et al. How- 
ever, the comparison of measured and calculated limb darkening is not a sensitive 
indicator of the photospheric absorption coefficient since the calculated limb 
darkening depends also on the photospheric model adopted. ‘This is clearly 
illustrated by the difference between the calculated limb darkening of Neven and 

- De Jager and Peyturaux (see Fig. 29 of Peyturaux) both using the same absorption 
coefficients but different models. Chalonge (19g), using his measured relative 
central intensities at different wavelengths in the visible and those of Mulders 
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(21) in the near infra-red and the limb darkening measurements in these regions 
of Abbot and his co-workers (22), concluded that an overall better fit with theory 
would be obtained if the free-free absorption was less by 25 per cent than given by 
(x). Our observations are not inconsistent with Chalonge’s (19) inference and 
the recent calculations of the free component of the absorption coefficient of H~ (8). 
However, they also imply a slower increase of absorption coefficient with wave- 
length in the far infra-red than k <A? and they do require significant modification 
of current photospheric models. 

Taking the bound-free (1) and the free-free (8) absorption coefficient of H~ in 
the near infra-red and the measured value of /(() (3) at 1-009 to be correct 
(the limb darkening measurements of Peyturaux here being used) it appears, from 
consideration of the spectral variation of the optical depth (Fig. 8) that Peyturaux’s 
values of (©) between 1-6 and 2-3 are systematically too low by approximately 
12 per cent while those of Mulders between 1-75 and 2-ou are too high and above 
2°2u are too low by approximately 12 per cent. 


0 us 


Fic. 8.—Experimental variation of the relative optical depth in the near infra-red at 
6000°K, compared with the theoretical variation of the absorption coefficient in this region. 
— — —Peyturaux's (3) values 
t+ + Mulder’s (21) values as published by Peyturaux (3) 
~ Theoretical curve of the relative absorption coefficient from bound-free (1) 
and free-free (8) at a temperature of 5500°K which is a good representation of the upper photo- 
sphere in this region. 
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ON TIDAL TORQUE AND ECCENTRICITY 
OF A SATELLITE’S ORBIT 


Gordon W. Groves 
(Received 1960 August 9)* 


Summary 

It has long been known, by a simple and general argument, that tidal 
dissipation within the Earth must slow its rotation, slow the Moon’s orbital 
motion and increase the Moon’s mean distance. A hypothesis is made 
whereby it can be shown, under almost as general conditions, that tidal 
dissipation must also increase the eccentricity of the Moon’s orbit. The 
exact rate of increase is not determined. However, it is shown that if the 
Moon were at one time quite close to the Earth the eccentricity would have 
been at most about 2 x 10~*; i.e., the orbit would have been nearly circular. 
This agrees nicely with Darwin’s idea that the Moon originated near the 
Earth. At the stage where, according to Darwin, the Earth and Moon each 
present the same face toward each other and revolve with a period of 55 
present days, the eccentricity will be about 0-2 to 0-4. 


Variation of an osculating element.—If « designates any osculating element 
of the Moon’s motion, its rate of change can be expressed by 
Here a(x, y, z, %, y, #) is determined, at any instant, by the position r= (x, y, 2) 
and velocity v = (x, y, #) of the Moon relative to the Earth at that instant. R is 
the disturbing function. It is convenient to define the operator 


a 


and to write the perturbation equation in the form 
x=Va.UR, (1) 
It is seen that Va is a vector in space whose value at the Moon’s instantaneous 
position indicates the direction in which a small impulse of given intensity will 
produce the maximum increase in «. As an example let us take « to be the total 
energy of the Earth-Moon system, 
_ 1_mgm 17 ymgm 
E te” (2) 
The three terms are the kinetic energy associated with the orbital motion, the 
kinetic energy of the Earth’s rotation and the potential energy, respectively. 
Here mand m, are the masses of Moon and Earth, respectively, y is the gravitational 
constant, w the Earth’s angular velocity, and J the Earth’s moment of inertia 
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about its axis, taken to be constant. The potential and kinetic energy of tidal 

motions need not be included in (2) since they are negligible compared to 

the other terms. The problem under consideration involves a change in E by 

an appreciable fraction of itself during the evolution of the Earth-Moon system. 
The total angular momentum of the Earth—-Moon system is 


H = lw (3) 


where h=r xv. All solar and planetary perturbations acting on the Earth-Moon 
system are neglected. In unperturbed motion h, and consequently w also, 
are constant. In perturbed motion h and w must vary in such a manner that H 
remains constant. Eliminating w between (2) and gn gives 


Rat 4 2" Aer xv 
2mM,+m Mo 
mym 
We then have 
= w xr). (4) 


The second term is greater than the first by a factor of about 30, and almost 
oppositely directed since the Moon revolves in the same sense as the Earth’s 
rotation. Since tidal dissipation converts energy into heat, the disturbing 
force cannot have a component in the direction of VE. Then in order that 
VE.VR<o the disturbing force must act in such a way as to increase the length 
of day, length of month and mean distance. Essentially the same argument 
is brought forth by Jeffreys (1959, p. 232). 


Eccentricity and semi-major axis.—Now let us take « to be the square of the 
osculating eccentricity. This can be calculated at any instant according to 


h2 
(5) 
where a is the osculating semi-major axis, equal at any instant to 
2\—1 
€ ) (6) 
ro 


and yp is written for y(m,+m). From (5) and (6) we have 


ve = (7) 


= (8) 
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Expression (7) is shown graphically in Fig. 1 for the case e=o-5. The rates of 
change can be calculated by (1): 


ie (:)"]v-ve + (9) 


(10) 


The disturbing acceleration VR is related to the tidal torque on the Earth L 
according to 


L=-—“ sx UR 
my +m 
from which we obtain 
verxL= ry. VR}. (11) 
Using (11) to eliminate v.VR from (g) and (10) gives 


de® _ 2a (2) (12) 


at pr | mom 


mgm 


2 (13) 


Scale of Vein 
units of 


Fic. 1.—The best way to push the Moon in order to increase the eccentricity. The vector 
Ve is shown at points along the orbit for the case e=0°5. 


Up till now no stringent hypotheses have been made. If the exact form of 
the disturbing force were postulated, the rates of change of a and e could be 
calculated from (12) and (13). If we abandon the attempt to determine the 
time rates of change and try only to evaluate the rates of change with respect to 
each other, a simple result can be obtained from quite general considerations. 
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Hypothesis I.—We shall assume that the in-line component of the disturbing 
acceleration, r,.VR, has a negligible effect on both a and e as compared to the 
perpendicular component associated with the tidal torque L. This may not be 
as unrealistic as it may seem at first glance. We have seen that the sense in 
which L must act is well determined from energy considerations; i.e., L.w =o. 
This condition must persist throughout all time that the Earth’s rotation is 
faster than the Moon’s revolution. On the other hand there is no general guide 
as to the sense of r,.V R, and the most natural thing to assume is that sometimes 
it acts in one sense and at other times in the opposite sense. The slowing of 
the Earth’s rotation by the accumulating action of the tidal torque has been 
observed. ‘There are as yet no observed phenomena attributable to the in-line 
component of the disturbing force. In any case we are not obliged to consider 
all effects simultaneously and so there is no harm in confining our attention to 
the effect of the tidal torque, as many previous discussions have done. 

Then dividing (12) by (13) gives 

2 2 
If the orbit is circular then expression (14) is always zero, and eccentricity is 
not spontaneously generated. If there is already some eccentricity, then the 
sign of (14) alternates, being positive at apogee and negative at perigee. To find 
the net effect it is necessary to make the following hypothesis. 


Hypothesis 11.—The magnitude of the tidal torque at any point in the orbit 
is assumed to vary only with the Moon’s distance. Let us take the plane of the 
orbit as the xy plane. Then 

v.rxL=-h.L=—-CN (15) 
where C and N are the z-components of h and L, respectively. In accordance 
with the second hypothesis we take 


N=] (Z). (16) 


Expression (16) does not represent the variation of N in detail, but gives, let us 
say, an average over many months for each particular value of r. Np is allowed 
to vary slowly with time (or semi-major axis) to take into account the varying 
constitution of the Earth over geologic ages. The greatest torque occurs at 
perigee, and the least at apogee, so that the function F will decrease with 
increasing argument. 

There is one important effect that this hypothesis does not take into account: 
the age of the tide; i.e., the time lag between the occurrence of perigee or apogee 
(or syzygies) and the corresponding maximum or minimum tidal range. The 
age of the tide depends on the relative phase lag for different frequencies. In 
the ocean, the age of the tide varies a great deal and is actually negative in some 
parts. A representative value, however, might be of the order of two days. 
Since this is short compared to the month we should not expect a large error to 
arise from our disregard of the dependence of tidal torque on the Moon’s 
longitude relative to perigee as well as on the distance. 

For an orbit which is not too eccentric, the quantity 7/a remains near unity 
and we can expand the function F in a Taylor series about this value: 


F(Z) = F(1) + ( 1) + 


aver 

at 

| 
| 
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Neglecting the subsequent terms, equations (12) and (13) can now be written as 


+ (-e)(2) - 


These expressions will now be averaged over a monthly cycle, during which a 
and e can be considered constant. The mean value of (a/r) raised to a power 
can be calculated from the well-known expansions for an elliptic orbit (see, for 
example, Smart 1953, p. 41): 


((2)’) = 1+ }n(n—1)e + O(e). 


The eccentricity is taken to be a small number. The first term of (17) and (18) 
in the power series in e is 
de? _ 2CNomyt+m 
= 
_ 2CNym)+m 


dt 


(19) 
Equation (19) can be integrated immediately to give 


where a, and é, are any pair of simultaneous values of a@ and e. 

It is of interest to estimate the numerical value of the exponent A in (20). 
Let us consider two cases which may give a somewhat unusually high and low 
value. First, if the tides on the Earth are linear, then the amplitude of the tidal 
bulge is proportional to r~* while the ee s torque on a given bulge goes also 
like r~*. This gives F~ (r/a)-* and A= Second, in the case of non-linear 
dissipation of energy by tidal currents in = ocean the drag is nearly proportional 
to the square of the velocity while the dissipation is nearly proportional to the 
cube. Ifthe current itself is nearly proportional to the inverse cube of the Moon’s 
distance, then this ‘‘quasi-linear’’ model dissipates energy ae r~*, Since 
torque is proportional to dissipation, we have F ~ (r/a)~® and A= We might 
expect that bodily tidal friction obeys the linear law more closely xn the ocean 
tide obeys the quasi-linear law more closely, but this is quite speculative. 


Conclusions.—According to (20) a circular orbit is stable with respect to a 
small initial eccentricity only if the tidal torque decreases with Moon’s distance 
no faster than the square root of distance. This is hardly likely in view of the 
two torque models discussed above, for both of which A is a large positive number, 
We should expect, then, that any small departure from a circular orbit will grow. 


36 
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If a, and é, are taken to be the present values for the Earth-Moon system, the 
following table of values can be constructed from equation (20): 


linear model 


quasi-linear 


model 
a (km) a/ay A=11/4 A=17/4 
Initial 10,000 e=2x10~* e=107° 
Present 380,000 I €=0°055 €=0°055 
Final 600,000 1°6 e=o'2 e=0'4 


“‘Tnitial’’ refers to Darwin’s (1908 or 1911) idea that the Moon originated 
quite close to the Earth, if not actually having once been part of the Earth. The 
initial value of a was chosen to represent. the Moon soon after its separation 
or formation. ‘‘Final’’ refers to Darwin’s described state in which the Moon 
and Earth each present a constant face to each other and revolve with a period 
of 55 present days. By either model the initial orbit would have departed 
very little from a pure circle. This tends to substantiate Darwin’s theory of 
the origin of the Moon. The detailed time history of the eccentricity as given 
here differs from Darwin, who found that even the sign of the change in 
eccentricity was governed by quite intricate conditions. ‘The exact manner in 
which Darwin’s viscous-earth tide differs from the conditions hypothesized 
here is not clear (to the author). 

The present rate of change of e can be estimated from the observed linear 
discrepancy and Kepler’s laws. The rate of change of a is estimated (see Munk 
and MacDonald 1960, Chapter 11, or Groves 1960) to be 10~‘cmsec~! or 
10 per cent of itself per billion years. Thus e increases at about 30 per cent 
to 40 per cent per billion years. Urey et al. (1959) have proposed a mechanism 
whereby tidal dissipation within the Moon’s body decreases the eccentricity 
and suggest that the Moon’s orbit was more eccentric in the past. It would 
be of interest to evaluate the magnitude of this effect in comparison to the opposing 
effect described here. 
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